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I. INTRODUCTION 



At high temperature T, the properties of the quark-gluon plasma are described in the 
context of massless QCD with a small (running) coupling constant which implies there is 



a hierarchy of scales g'^T, n a nonnegative integer_ 
determined to lowest order in the coupling P, [T 



The quasiparticle spectra have been 
. There are transverse gluons and 



longitudinal ones (plasmons), ordinary quarks and plasminos In ^ was recognized the 
problem of gauge dependence of naive one-loop-order dispersion relations for slow-moving 
particles. A systematic reorganization of the perturbation became necessary and a proposal 



was made in 
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IJ]. It consists in summing the so-called hard thermal 



loops (HTL) into dressed propagators and vertices [if 

Hard thermal loops form a lowest-order approximation of massless QCD at high temper- 
ature describing slow-mo ving qu asiparticles. They have been encapsulated in a generating- 

and once a connection to the eikonal of a Chern- 



moving quasipartic 



functional formalism 
Simons gauge theory was made 



HQQ, a hydrodynamic approach 12^,12^12^ showed 
that they are essentially classical. This suggests that 'true' quantum effects would effec- 
tively arise at next-to-leading order in HTL-summed perturbation, and one question to ask 
is whether this perturbation is consistent and reliable. More particularly, is a next-order 
term really a correction or of the same magnitude as the one considered? Also, does the 
HTL-summed perturbation take account of and describe all p ossible effects of the theory? 
Indeed, note that the functionals in the literature 24 , 2^ bfil ] generate only the hard ther- 



mal loops and not the full quantum field theory organized in a HTL-summed perturbation. 
Finally, are actual calculations safe from potential infinities? 

With this in perspective, many quantities have been examined in next-to-leading order 
calculations. In gene ral , finite and consistent results are obtained when fast-moving particles 



are involved 
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3G 



31 



Q, 33, QQ, 



but difficulties arise when slow- moving 



quasiparticles are dealt with. For example, the determination of the production rate of soft 
real photons from a hot quark-gluon plasma 0] has revealed that the hard thermal 
loops introduce non-cancellable coUinear divergences when external momenta are light-like. 
It is shown in another approach ^| that it is not even the hard thermal loops that dom- 
inate at leading order. The problem of collinear divergences is looked into in 38] and an 
improved gauge-invariant effective action is proposed that takes into account asymptotic 
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thermal masses. Further work 



39|,|4C, 



4l| on the production of photons from a hot quark- 



gluon plasma using this improved action shows that (i) two and three-loop graphs are of 
the same magnitude as the lowest order term; (ii) though infrared divergences associated 
with ultrasoft transverse gluons cancel in graphs with what is termed 'abelian topography', 
due to a sort of a finite-temperature Kinoshita-Lee-Nauenberg theorem, the result depends 
strongly on the introduced magnetic mass. 

Regarding the infrared sector more specifically, there are strong indications that here too, 
problems may arise when pushing the HTL-summed perturbation to next-to-leading order. 
Indeed, it is well known that at lowest order gT, static chromoelectric fields get screened 
but not the static chromomagnetic ones. These are believed to get so in nonabelian gauge 



theories at next-to- leading order g'^T 
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44| . the so-called magnetic scale. It is also 



believed that the determination of the chromomagnetic screening length, the inverse of the 



magnetic mass /i, is nonperturbative 
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47| . Therefore, it is quite possible that non-screened 



static chromomagnetic fields spoil some, if not all, next-to-leading-order calculations. For 



example, the works 



4a, |42| discuss next-to-leading-order nonabelian Debye screening in 



HTL-summed perturbation and it is explicitly shown that the results depend strongly on 
the type of infrared regularization used. 

There also seems to be difficulties in the infrared when dealing with next-to-leading-order 
dynamic quantities. The first non-trivial such quantities to investigate are the slow-moving 
quasiparticle damping rates. Indeed, though they come from one-loop calculations with 
fully dressed propagators and vertices, they only involve the imaginary part of the one-loop 
dressed self-energies, something that simplifies matters considerably. The first work to cite 
is j3] which determines the damping rate 7t(0) for non- moving transverse gluons and finds 
the finite and positive result: 

^tiO) = 0.088N,gX (1.1) 
where A^^. is the number of colors. Another result is that of the damping rates 7+(0) for 

n n 

non-moving quarks found independently in j49|| and j50|| which read: 



7±(0)=aoC//T, 



:i.2) 



where Cf = {N^ — 1) /2A''c and ao is a finite positive numerical constant that depends on Nc 
and the number of fiavors Nj. However, in contrast with these finite results, the damping 
rates for on-shell slow-moving 



the presence of a logarithmic divergence in g: 

7(p) ag'^Tv{p) \ng, 



(1.3) 



where a is a positive numerical constant, v{p) the quasiparticle group velocity and p = |p|. 
Note that if one takes the limit p — > 0, one cannot recover neither the finite result ()1.1|) nor 



fll.2|) . Another estimate is done in 53] and a result similar to ()1.3p is found: 



7(p) ^-^^v{p)\ng. (1.4) 

The infrared problem in the damping rates is more emphasized in |^ , where a direct 
and explicit calculation of the damping rate 7/(0) for non-moving longitudinal gluons in the 
sole framework of the HTL-summed perturbation shows that this latter is infrared sensitive: 



7/(0) = ( 7^;rT^ - 2.471908492 ) (1.5) 



250 \ g^N^T 

— — - 2.471908492 ^ 
277r2/i2 J 247r 

where ft is the magnetic mass /i, taken as an infrared regulator, in units of the electric 

thermal mass. This peculiar result contrasts with the expectation that at zero momentum, 

there must be no difference between longitudinal and transverse gluons, and so the two 



so, 



corresponding damping rates ought to be equal |48|, a statement emphasized in j55|. A 
another direct and explicit calculation of the damping rates 7±(p) for ultrasoft | 
quarks suggests that these too bear infrared divergences. 

The purpose of the present article is to determine the damping rate 'jt{p) for a transverse 
gluon with an ultrasoft momentum p. The aim behind is an additional attempt to un- 
derstand the infrared behavior of massless Quantum Chromo dynamics at high temperature. 
The calculation is done in the strict context of next-to-leading order HTL-summed perturba 



tion theory as developed in 



m 



IJ] . In the back of our mind is the persistent 



question whether this scheme is reliable and efficient enough when it comes to describe the 
properties of the quark-gluon plasma, most importantly beyond lowest order. We will seek 
an expression for '~ft{p) in an expansion in powers of the ultrasoft external momentum p. We 
recover the first, finite and positive, coefficient p.ll) found in All coefficients of odd 
powers of p vanish. We determine explicitly the coefficient of p"^ and find it sensitive in the 
infrared, see ()4.43p . This finding relies on the way the calculation is performed and much 
of the discussion we carry orbits around this point. This may be felt slightly unfortunate a 
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situation, but it should be seen as an additional reflection of the degree of difficulty to carry 
analytic estimates of physical quantities beyond lowest order directly in QCD. 

In the next section, we give the definition of the damping rate in the context of HTL- 
summed perturbation. The problem amounts to calculating the imaginary part of the 
transverse-gluon one-loop self-energy with fully HTL-dressed propagators and vertices. We 
carry this calculation in section three. The analytic results of this section have already been 
reported in j58i] and so, this part of the work is not new; we merely give here substantially 
more details regarding the intermediary steps. The evaluation of the analytic integrals is 
carried out in section four. It necessitates the disentanglement of the infrared- sensitive pieces 
from the finite ones, quite delicate a task. We choose the magnetic mass fi as the infrared 
regulator. It also requires the handling of potential divergences coming from light-like mo- 
menta. We explicitly show that these latter can be circumvented without having recourse to 
the improved HTL-summation of js^, whereas the infrared divergences persist. All of this 
section is new. 

Section five is devoted to a detailed critical discussion of the result we obtain and the 
method we use. First we explain why it is admissible to regulate the infrared sector using 
the magnetic mass fi, and assuming that all of massless QCD is in the HTL summation 
scheme, we impose the equality between 7i(0) and 7t(0) to obtain a rough estimate of p,. 
Various orders of magnitude are discussed and we come to the conclusion that the HTL- 
summed perturbation may be useful well into the quark-gluon- plasma phase. Next we 
address the issue of the equality that we do not obtain in between the damping rates 
for the longitudinal and transverse non-moving gluons, in view of the work in which 
a Slavnov- Taylor identity for the gluon polarization tensor in Coulomb gauge is derived 
and, when formally applied to the next-to-leading order self-energy, equality is found. We 
show that this identity, though important in itself, is not of direct help to us. We then 
compare our result with other estimates in the literature, 0, Q , and argue that there is no 
contradiction between what we find and these works: the regions of the external momentum 
p are not the same. To indicate that our method of calculation has little to do with the 
occurrence of infrared sensitive coefficients in ultrasoft physical quantities, we carry a similar 
calculation in the context of scalar quantum electrodynamics where it is possible to avoid 
the early expansion in powers of the external ultrasoft momentum. We show that results 
with early expansion and without are the same. We conclude the discussion with few final 
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comments. 

One appendix contains the small-momentum expansion of the gluon on-shell energies, 
the residue and cut functions involved in the spectral densities of the dressed propagators, 
together with their first and second derivatives. We should mention that only the first few 
terms of the expansions are displayed; in actual use, many additional terms are needed. A 
second appendix gives the results of the terms contributing to the imaginary part of the 
transverse-gluon one- loop HTL-dressed self-energy we have not detailed in the main text. 



II. GLUON DAMPING RATE IN HTL-SUMMED PERTURBATION 



The calculation is performed in the imaginary-time formalism where the euclidean mo- 
mentum of the gluon is = {po, p) with = Po + and the bosonic Matsubara 
frequency po = 27rnT, n an integer. In this formalism, real-time quantities are obtained 
via the analytic continuation p^ = —iuj + 0"'', where uj is the energy of the gluon. A momen- 
tum is said to be soft if both uj and p are of order gT] it is said to be hard if one is or both 
are of order T. The three-momentum p of the on-shell gluon is said to be ultrasoft if p is 
much smaller than gT, say of the order g'^T and smaller. The calculation is carried in the 
strict Coulomb gauge = where the separation between the longitudinal and transverse 
cornponents is straightforward. The HTL results we quote in this section can all be found 



m 



)onents is straigntiorward. 

mum 



In the strict Coulomb gauge, the dressed gluon propagator *A^,^ (P) has a simple struc- 
ture. It is given by *Aoo (P) = ^A^ (P), *Ao. (P) = and *Ai, (P) = {5^, - p^Pj) *At (P) 
with: 

where 511; (P) and SIlt{P) are hard thermal loops given by: 

(2.2) 



OTKP) = SmlQ, ( ^ j ; 5n,(P) = ^m; 



The Qn are Legendre functions of the second kind and nig = ^/N\^'^^TjJJ2 gT / ?> to low- 
est order is the gluon thermal mass. The poles of *At(i)(— icj, p) determine the transverse 
(longitudinal) gluon dispersion relation. These write: 

, ut + k 2k(-2euj^ + 2k^ + 3mlu^) , ui + k 2fc (Sm^ + 

u = ^2 ' 1^ u = ^ • (2.3) 

uJt — k Srrig Ut [oof — /e ) ooi — k Srrig lji 
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At this lowest order gT in the dispersion relations, the on-shell gluon energies are real and 
no damping occurs. To get the damping rates to their lowest order g^T in HTL-summed 
perturbation, one has to include in the dispersion relations the contribution from the next- 
to- leading-order self-energy *n which has a more complicated structure than that of the 
corresponding hard thermal loop. It satisfies the less restrictive identity P^P^, *n^j,(P) = 0, 
which means *n depends on three scalar functions. The advantage of the strict Coulomb 
gauge is that only two of these are relevant, namely: 



The transverse-gluon dispersion relation reads: 



(2.4) 



(2.5) 



The transverse-gluon damping rate is defined by '~ft{p) = —l^aQt (p)- It is (yf-times smaller 
than ujt{p), and so we have to lowest order: 

lm*Ilt{-iu!,p) 



lt{p) 



2.6) 



2uj + d^6Ilt{-iuj , p) 

where d^^ stands for d/du. Using the expression of 5Iit in ()2.2|) . we have for soft momenta: 

2 



lt{p) 



10 Vm, 



Im*nf(— iti;,p) 



(2.7) 



uj=u)t{p)+iO+ 

Hence, determining 7t(p) amounts to calculating the imaginary part of the next-to- leading- 
order transverse-gluon self-energy. 

In HTL-summed perturbation, the next-to-leading-order gluon self-energy in the strict 
Coulomb gauge is given by j3 ]: 



m^-{P) = -^Trsoft [T^^"'^(P, -P, K, -K) *^x.{K) 

+ *r'^A.A(_g^ *AAA'(i^) T^''^'^'(-^, P. -Q) *^.'.{Q) 



(2.8) 



K is the loop momentum, Q = P — K and Tr = J - — — with kg = 2mTT. The subscript 

ko (271") 

'soft' means that only soft values of K are allowed in the integrals, which implies that both 
propagators and vertices must be dressed; this is already indicated in ()2.8p . The dressed 
vertices *r are of the form: 



(2.9) 



where T^^^^ is the tree three (four)-gluon vertex and ST^^^^ is the corresponding hard thermal 
loop. We have: 

with 5* = (i, s) and Qg the solid angle of s. 

To get an expression for the imaginary part of *Ilt, the next-to-leading-order transverse- 
gluon self-energy, we have first to perform the above angular integrals, then perform the 
Matsubara sum, then analytically continue to real gluon energies po = —iuj + 0"^ and imple- 
ment the on-shell condition. Last is to find a way to perform the integration over the gluon 
loop three- momentum k. However, it is technically difficult to follow the above sequence 
of operations. Especially the integration over Vt^ in (j2.1Uj) and (|2.1ip . To render this latter 
feasible, we expand *n in powers of the external gluon momentum p/mg, a small parameter 
when p is ultrasoft. The Matsubara sum is performed using the spectral representations of 
the different quantities involved. The subsequent integration over Qj. is analytically done 
and the remaining integrals have to be performed numerically, once infrared-sensitive pieces, 
if any, are extracted. The rationale behind this procedure is as follows. Expecting infrared 
sensitivity, we introduce from the outset an infrared regulator which remains fixed for the 
rest of the calculation. It can be viewed as physically representing the magnetic scale g^T, 
or better the magnetic mass /i. The integration J^^°° dk is therefore replaced by /^°° dk. 
This means that the loop momentum k is always greater than fi. As we said, we always 
regard in this work p as ultrasoft, i.e., p < fi < k. This condition allows the expansion of 
quantities function of g = |p — k| in powers of p, quantities like 1/QS and *A{Q). We have 
already discussed the range of validity of this expansion in 0] and will come back to this 
issue later in section five. 

III. ONE-LOOP HTL-SUMMED GLUON SELF-ENERGY 

From now on, we set rrig = 1. This will simplify writing the intermediary steps as well as 
the expressions we obtain. From ()2.8|) and the structure of the gluon propagator *A^j, given 
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\P, -P, K, -K) 5mn - kmkn * MK) 



in the text before ()2.ip . we have the following explicit expression for *Ilt: 

^ fco (27r) 

+ *r°*°(-Q, P, -K) *MK) T°^°(-ir, P, -Q) *Ai{Q) 
+ *rO.™(_g^ (^^^^ _ *Y"j%-K, P, -Q) *AKg) 

+ *r™'-(-g, p, -ir) (drs - kk) *Mk) *r^^{-K, p, -q) (5„^ - *Ai(g) . (3.1) 



There are six contributions: two from the four-gluon (4(7) vertex and four from the three- 
gluon {{"ig) vertex. Each contribution has to be calculated separately. We will show the 
details of how we get the Sf/ZZ-contribution where the two propagators are both longitudinal. 
We will comment afterwards on how to work out the other contributions. 

A. Solid-angle integrals 



Using the expression of the dressed S^f-vertex, we have: 

•n.3,„(/') = ^rj; / ^ [(p - 2k)^ - (p - 2k).p' 



+ 6 
+ 9 



dVLs (p — 2k). s — (p — 2k).ps.p / ik^ iqo 

17 PS " QS 

dQ.si f d^ls2 S1.S2 - Si-Ph-P f iko iqo 



ikn 



*Ai{K)*Ai{Q). (3.2) 



47r J Att PS1PS2 \KSi QSiJ \KS2 QS2, 
We first work out the contribution that involves no integration over Qg- We have (p — 2k)^ — 
(p — 2k). p^ = 4fc^ sin^ -ip with = ( p, k) . To integrate over the solid angle fi^, we expand: 



,0 f sin^ ip ^ 2 , o2 

l-pcos'tpdk + Y 1 — — "^fc + cos V Ok 



^A,(,)(go,fc), (3.3) 



where dk stands for d/dk. The solid-angle integrals become straightforward and we find: 
d^k 



— [(p - 2k)^ - (p - 2k).p^] *A,{K) *A,{Q) 

KQ 



k^dk*Ai{K) 



1 + 1(^4 + 1^1+... 



'Ai{qo,k). (3.4) 
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Next we work out the term involving one single integral over Qs- It is sufficient to 
concentrate on the piece with ik^/KS because the other one, that involving iqo/QS, is in 
fact equal to the first one. To get a manageable expression for the solid- angle integral, we 
use the following expansion: 



1 



1 



p.s p.s 



Wo 



Also, it is best here to measure the solid angle Vts 
9 = fk, . We can then write: 

iqo 



+ ... . (3.5) 
, ip) with respect to k such that 



dVls (p — 2k). § — (p — 2k).ps.p / iki 



-Ak sin ip 



PS \KS QS_ 

iko f dQg sin ip cos 6 + cos ip sin 6 sin (p 



(3.6) 



ipo J 4:71 iko + k cos 9 [ ipo 

The angular integrals are performed using the relation p.s = cos ^/^ cos — sin?/'sin0sin(y9. 
We get: 



d^s (p - 2k). s - (p - 2k).ps.p / iko_ _ iqo \ _ /-^ _ 2\ 

Air PS \KS QS) ~ ipo^ ^ ^ 



2(1 — —^Qok 



. P 
H X 

iPo 



iko 



^(l-3x^ 



1 + 3|) Qo. 



p_ 

pI 



(l-2a:^) + (l-5x^)| 



(1 - 5x^) ^ ) Qo. 



+ . . . 



(3.7) 



where x = cosip and Qok = Qoijr)- We put this expression back inside the integral over 
d^k and use ()3.3j) to perform the integral over Qk- We get: 



ko 



ko 



Arc 

+ 00 



PS 



KS QS 



k''dk*AiiK)'-^ 
n Wo 



iko. 
k 



1 + '^'n 

+ —j^^Ok 



+?(3^.(i.3|)gJa. + p^(i-^^(|a..ia| 



+ ... 



*Ai{qo,k). (3.8) 



Next we look at the terms with two solid-angle integrals. They are actually equal to 

koqo 



-18TJ2 

ko 



— *Ai(K)*Ai(Q) [ /'c^^]s2 S1.S2 - si.pss.p 
Air I Att 



Att 



PSi PS2 



KSi KS2 KSi QS2 



and we ffist look at the piece containing ko/KSi KS2', the other one will require more labor 
as we will see. The additional difficulty here is the coupling between the two solid-angle 
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integrals introduced by the dot product S1.S2. In order to be able to carry forward, it is best 
to write S1.S2 — Si.ps2.p explicitly in trigonometric terms, which will automatically decouple 
the double solid-angle integral into the following form: 



4:71 



+ COS^ ijj 



dVls^ S1.S2 - Si.pS2.p 



47r 

dfls sin 6 sin 



An PSKS 



PSi PS2 KSi KS2 
2 

+ 2 cos ip sin ip 



dQs cos 6 
Att PSKS 

dQg sin 6 sin ip 



sin^ -ip 

dQs cos 6 
An PSKS I An PSKS 



(3.9) 



The problem reduces therefore to carrying out only a single solid-angle integral each time, 
which is straightforward using ()3.5|1 . We have the following intermediary results: 

dQ^ cos6 1 



An PSKS ipok 



ikQ px iko 

-L ^vofc + -. r- 

k tpo k 

F 2k 



1 

i '^^Ok 



o 2 1.2 



X 



+ (1-3.')|)q 



+ 



(3.10) 



dQs sin 6 sin ip sin ip 



An PSKS 



ipok 



P ( iko 



2ipo 



^+ 1 + ^ g 



A;2 



iOk 



p^x 



3 A:2 



k 



1 + 



ftp 
A;2 



Ofc 



+ 



(3.11) 



whereas we have / 



47r PSKS — ^ ■ Putting these back into ()3.9|) and organizing the prod- 
ucts and squares in powers of we get: 



j dQs, j 



dils2 h-h - Si.ps2.p 



/ 



An 



PSi PS2 KSi KS2 



p^k^ 



px 



k 



i + 3^)gofc 



ik 



iko 

i '^^Ok 



kn 



k 



^ + 1 + ^ Q, 



fc2 



iko. 



ka\ ikn 



+ (1 - f «o. j (i - (1 - 7r) g - (^1 - 4x= + (1 - 7.^) ^ j ^0„. jj + ..J (3.12) 

We put this expression back under J d^k and use ()3.3j) to perform the integral over the solid 
angle Qk- We get: 



ko 



d^k 

An 



Ai{K)*Ai{Q) 



dQ 



An 



dn 



S2 S1.S2 - Si.pS2.p 



An 



PSi PS2 KSi KS2 



+00 



+ 1 



ikn 



X 3- 



iko 



Ok 



1 + 3 



dk*Ai{K) 
) _ 2-0 - 



1 

-1- — '^^^^ 

A;2 



5p2 



iko 



+ 1 + 



Ofc 



iko 



1 — —^Qok 



dk + dl 



k 
+ 



'Ai{qo,k). (3.13) 
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Note that generally, odd powers of p do not cancel in the intermediary steps until the 
integration over Qk is performed. 

The piece in the double solid-angle integral in (j3.2|) containing k^q^/ KS1QS2 deserves 
more attention. Because of the presence of both k and q, it is most suitable to measure the 
solid angles with respect to p. Then we have: 



An 



S1.S2 - Si.psa.p koqo 



An PSi PS2 



+ 



KS,QS2 
dQs sin 6q sin 



dflg sin ^0 cos ipo f dflg sin 6q cos (Pq 



An PSKS 
dQs sin ^0 sin po 



An PSQS 

(3.14) 



An PSKS J An PSQS 
where we have added the subscript zero to indicates that the solid angle Qs = {^o, V^o) is 
measured with respect to p. Each solid-angle integral can be carried out separately using 
the following identities: 

cos ^0 = COS COS 6 — sin i/j sin 6 sin ip ; 



sin 6q sin po = sin ip cos 6 + cos sin 6 sin cp ; 



sin 9o cos po = sin 9 cos p . 

We get: 

/dfls sin 6*0 cos po 
An PSKS ^ ' 

using inrTHll and the fact that / ^ %Ik7 = 0- ^he other integral in l/KS is: 
dQs sin ^0 sin po sin ip 



(3.15) 



(3.16) 



An PSKS 

p2 



iko px 
tpok I k ZtpQ 

^(l-2x^) + (l-5x^)M_^_3x^ 



3^+11 



3^ 1 Qok 



(3.17) 



We have similar results for the integrals with 1/QS in ()3.14|) : we need only to replace K by 
Q and the angle ip = (p, k) by x = (p,q). We then multiply the two expressions and put 
back the result under J d^k/An. But here, the integral over Qj^ is not straightforward yet 
because the integrand still depends on q and the angle (p, q). Hence a further expansion is 
necessary, but instead of expanding Qog directly, it is most suitable as we will explain why 
shortly to express first Qofc and Qoq in terms of *Aj'^{K) and *A^^{Q) respectively, using 
andQ. We have: 

(3.18) 



t A -1 



e + 3 1 



iko 
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and a similar result when replacing K by Q. We have written *Aik for *A;(i^) for short. 
The reason behind this preliminary replacement is to eliminate the occurrence of products 
of more than two functions necessitating a spectral decomposition. Such undesirable prod- 
ucts complicate unnecessarily the subsequent extraction of the imaginary part. After the 
replacement (|3.18|) is done and apparent simplifications made, we obtain: 



fco 

sin ip sin x 



fdVts^ /"c/fi^^ S1.S2 - si.pss.p koQo 
47 ^''^ / ^ 



47r PSi PS2 KSi QS2 



fco 



d\ 

47r 



*A 



IK 



X 



-koqokq 



pxqo 
Po 



{3 + e + 3kl)q-f^ 



-2^ (1 - 5x' 



k 



+ (1 - 3x^) k^+{l- 5x^) kl)q + ^{3 + e + Sk^) (3 + + 3g^) 



+ 



*A^Q, (3.19) 



where y = cos x- We only need now a final expansion of functions of q in powers of p. We 
use the additional relations: 



qcosx = P — kcosip; qsinx = —ksinip , 



(3.20) 



and the following two useful expansions: 
, = ,^l_,| + ±(l_,^)|! + ...); J_ = J_(i + 2.|-(l-4.^)|; + ...j. (3.21) 

Inserting all this back into ()3.19|) and performing the integral over dQk which becomes 
feasible now, we obtain: 



,Qrr^f^^^*/K *A /"c^^si /"c^^sa S1.S2 - si.psa.p koqo ^ v-^ /■+~ P 



X 



-koqok^ + 



5Po 



2fcogofc' + - (3 + A:' + ?,kl) (3 + + 3ql) 



+qoPok {3 + k^ + Sk^ - 2kopo) dk - -koqoplk^d] 



'MQo,k). 



(3.22) 



B. Sum over Matsubara 



Everything is ready to perform the sum over the Matsubara frequency kg. It is worth 
mentioning that the main reason we went through the above steps is to get iko (and iqo) 
appear only in the numerator of fractions. This way, to perform the sum, we only need to 
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use the spectral representations 



59 



60| of the following quantities: 



Qo{iko/k) = -- 



l/T 



l/T 

dre 



+00 



ikoT 



dujpt,i{k,uj){l + n{uj))e '^^ ; 
dujQ{k - \uj\) {l+n{uj))e~ 



(3.23) 



where n{uj) = 1/ (e'^/"^ — l) is the Bose-Einstein distribution and pt^i{k,uj) the spectral den- 
sities given in f I4.2p below. After we replace *Aik and Qofc by their spectral representations, 
we perform the integrals over the imaginary times. As we mentioned in the text after ()3.18p . 
preventing the occurrence of the product of more than two functions necessitating a spec- 
tral decomposition will ensure the occurrence of at most two (imaginary) time integrals. 
One such integration yields a delta-function and the subsequent one an energy denominator. 
Thus, terms have been arranged in such a way to yield at most one energy denominator. 
Only now can we perform the analytic continuation ipQ ^t{p) + ^0^. But just before 
this, every e~ has to be replaced by 1 except in the energy denominators. At each time, 
we are left with two frequency integrals together with the one over k. Terms involving one 
integration over imaginary time are real and the imaginary part of those involving two is 
obtained using the known relation 1/ (x + iO^) = Pr (1/x) — i7r5(x) applied to the energy de- 
nominators; Pr stands for the principal part. This technique is applied to expressions (j3.4j) . 
()3.8|) . ()3.13|) and ()3.22|) . We perform the sum and obtain after few steps of straightforward 
algebra: 



247r 



dk 



duo 



duo' 



18A; 



Quo^ 



uouo 



+ 



p 



2 r 



3k' 
2uouo 



3k f 



uo 



kuo 
2 



^k^ I ,0 0\r^/ ( ^ ^'^\ ^ rs , 

+ — (13 + 3fc2 + lOcu - 9cu') pidkp'i + 6- 1 + 2— - 3— Qdkp\ 
uouo V k J 



54fc 



2^4 o 2 

+ — (2 +5C.) pidlp\ +—^dlp[ -^-^pip'A 

uouo 



t^pA 



S. (3.24) 



ku' uouo"^"^'^ kuo' 

The notation is as follows. The quantity pi stands for pi{k,uo) and pj for pi{k,uo'). O is 
Q{k — \uj\) and 6 is 5(1 — uo — uo'). In the above relation, we have used the approximation 
n(u;*^'^) ^ T/uo'^''^ j3] because only soft values of uo^'"^ are to contribute, as dictated by the 
HTL-summed perturbation. Also, the last two terms in ()3.24|1 come from an expansion 
of 5 [uot (p) — uo — uo'] in powers of p. Finally, we note that many terms from intermediary 
steps drop in the end since they do not carry an imaginary part, i.e., they involve only one 
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imaginary-time integral as explained. The first two terms in ()3.24|) are those found in j48|. 
The terms in the contribution are new. 

What remains to do is to perform the integrals involved in (j3.24j) . How this is done is 
shown in the next section. Before this, we comment on the other contributions to Im*nt(P) 
and give the corresponding results. 

C. The other contributions 

The other contributions to Im*nt(P) are worked out along similar lines. Concerning the 
3(7-contributions, the intermediary steps of course are much longer. This is because of the 
presence of the transverse projector {dmn — kmkn^ in the three-gluon It and tl contributions 
in ()3.H) and the additional one {6rs — QrQs) in the three-gluon tt contribution. The trigonom- 
etry in the numerator becomes less simple and the expansion in powers of p quite intricate 
at times. For example, for the It contribution, to reduce terms involving a product of more 
than two functions requiring a spectral decomposition, we need to use, in addition to ()3.18j) . 
the relation: 

But in order to avoid dividing by a frequency (or go) while using this relation, something 
that would spoil the efficiency of the spectral-decomposition method, the Legendre function 
Qok has to appear in the different expressions in the form ^^Qok, which turns out to be not 
always the case. In contrast, Qofc in the 3gll contribution we just worked out had to always 
appear only in the form ^Qok, see ()3.18|) . This latter requirement is less stringent and 
it was always satisfied there, which meant there were no additional (technical) difficulties. 
What we use to remedy to the situation here is to group different terms together to obtain 
a cancellation of undesirable Qoq's. While for the tl and It contributions this is more or less 
straightforward, the cancellation is quite tricky for the 3gtt contribution since it necessitates 
grouping terms coming from the expansion of the zeroth and first orders when using ()3.5p . 
terms seemingly not connected together. These intermediary steps work, but only case by 
case: we do not have a general rule regarding this point. These steps are quite lengthy and 
cumbersome to report on with any degree of detail. We limit ourselves to giving the final 
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results, similar to ()3.24|) . We find for the tl and It contributions: 



A8tt 



dk 



-oo ^ 
2 



— oo 
2 



duj' 



UJ' 



-18 {e - u^Y p,p[ 



+ (8 + rSA;^) - (+140 + 78k'^ + 90k^) + (48 - 180P) uj 



- (31 - 126A;2) + 102cu^ - 54cu^] pjpj + 



2k 



3 

+ 3A; ) — -^ti; 



+ 



3 15 



k^ 



UJ 



57 21 
2k^~l^ 



24 

-6kuj + — cj^ 



18 

fc3 



a; 



Gp; + [-12/c-72P+(-24/c -12^;^ + 48P) uj + (-12A; + 144A;3)a; 

+ (12A; - 96k^) - 72ku^ + 48ku^] ptdup'i + [69A; + lAk^ - 3k^ 
- (54A; + 16k^ - 6P) u + {2k + 14P) uj^ + {8k - 12P) 

-llkuj^ + 6kuj^]pidkp't + 

edkp'i - [8A;2 + 16k^ + IGk'^u 



9 ^ 2 3 3 18 4 3 5 9 f 



+ 



9uj - 18uj' 



9 3 18 4 



+ (8A;2 - 32fc^) cu' + lefc'cu"] ptdlp'i + [SOk^ + Ak^ - 2k'^ - {24^ + 8k^) 
- {4k' - 4A;^) uj' + 8euj^ - 2^^^] pidU - 3-^ {k^ - oj'f Qdlp[ 

rh 



+ 6^ {e - uj') edlp\ + 54fc2 {k' - uj'fptp'id^ + 9^ {k' - uj'Y Qp'A 



- 18y {e - uj') Qp'id^ 
k 



+ ... 



6. 



(3.26) 



Note that the two contributions tl and It are equal as it should be. The above expression is 
rather long, and this is partly due to the fact that the functions corresponding to these two 
contributions do not benefit from a symmetry between uj and uj', present in the functions 
corresponding to the // and tt contributions. Indeed, though the algebra for the tt contribu- 
tion is by far the most tedious, the use of this symmetry, after the obtainment of expressions 
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like those of ()3.26|) above, renders the corresponding final result relatively simpler. It reads: 



247r 



dk — 



UJ 



'duj' 



UJ' 



36 (p + cucu')'ap;-6^ {e 



+ 



174 



-18A;2 + 122A;" - (l44 - 174A;2) ujuj' + - 138^ o^V^ 



'iuj 



+ 



32P + 20k^ + {112k + UAk^) cu + - 36k - 116k^^ cu^ 



108 



+ 232A;^ tu'' + 



12 



+ 156A; )cj^ + ^tu^-^cj 
J k k 



108 



60 



+ 12 (l + - (12 + 32A;2) tu^ - 12w^ + 12w^] pt^^pj 



UJ 



- 3^^ - uj') QdU - 108 (A;2 + uuj'f p,p',d^ 



k^ 



+ 18^ {k' - u') Qp'A 



+ 



6. 



(3.27) 



The two 45f-contributions require less labor but need few additional ingredients. Consider 
for example the four-gluon longitudinal contribution. Using result ()2.1ip . we have: 



kg 



d^k 
(2^ 



4-3 



dQs (1 - s.p2) f i{po- ko) i (po + ko) 



An PSKS \{P-K)S {P + K)Sj_ 



In addition to the expansion (j3.5j) . we use: 



1 



±1 



[P±K)S K±S 



p.s p.s^ 



K± = {ko±po,k) 



(3.28) 
(3.29) 



K±S K±S^ 

Putting back in (j3.28p . we encounter terms of the type 1/KS K^S"', with n = 1,2 and 3. 
They are handled using repeatedly the relation: 



±1 



KSK±S ipo 



1 



K^S KS 



(3.30) 



The rest carries straightforwardly until we obtain: 



47r2 



fco 



+ 00 



dkk^ 



iko 
k 



5Po 



E 



k, 



Oe 



2k 



Oe 



ik, 



ikn 



Kf poK^ plk 



Pok 



*A 



IK ■ 



(3.31) 
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To perform the Matsubara sum, we need the additional spectral decompositions: 

1 /-vr 



duj{l + n{uj)) e (lu) 6{uj^ - k^) e 



2\ -U}T 



Jo J -oo 

1 /-l/T /.+00 

^ -/o i-oo 



where e (u;) is the sign function. Few steps more and we obtain: 



247r 



dk / du 



dJ 



18k , 



(jj' uo' uo' 



+ ... 



(3.32) 



5. (3.33) 



The transverse ^-contribution carries in as much the same way. We obtain for it: 



Im*nt4gt(P) = _ ■ / ciA; / rfcj / dcu' ^ 



247r 



p2 

+ 5 



12 Icj 



kuo' 

, 6(4cj- 1)^ , 18(A;2 -cj2) ^ 



a; 



©pis'. 



+ 



5. 



(3.34) 



IV. THE DAMPING RATE FOR TRANSVERSE GLUONS 



The damping rate for ultrasoft transverse gluons to lowest order g'^T in HTL-summed per- 
turbation is given by ()2.7|1 where the imaginary part of the HTL-dressed one-loop transverse 
self-energy is given by: 



Im*ni(P)= im*n,3,.,(P) + 5^im*nMg.(P). 



(4.1) 



The different contributions are given in (jT^ . (jT^ . (HrTTj) . (pO!?|l and ( ICTl) . These 
expressions are in the form of triple integrals which ultimately are going to be performed 



numerically. The expressions o: 



propagators cirG clS follows 
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the dispersion functions corresponding to the dressed gluon 



pt^i{k,uj)=iti{k) [5{uj-uJt^i{k))-5{uj + uot^i{k))]+PtAk.^)Q{k~\^\), (4.2) 



with the residue functions given by: 



3cu2 - iuj^ - k^ 



UO (Cj2 - k^) 



UJ=LjJt(k) 



P (3 - cj2 + P) 



(4.3) 
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and the cut functions by: 

Pt{k,uj 



3a^ 

3u 



[k^ a;2 + 1^ (l + ^^xx. ^ T6fc6 



In 



k+u 



2k 



(3 + A:^-ilnS) + 



(4.4) 



Recall that the thermal gluon mass rUg is set to one. There are different types of terms 
involved in (HT^ . (HT^ . (ITTfll . (HUIHl) and (ESI)- We will show how we carry through 
using a generic term of each type. Because of lengthy expressions, we will use as much 
compact a notation as permissible in such a way that the meaning is always clear from the 
context. 



A. Integration 

The first type of integrals we have to deal with is the one that involves a pp contribution. 
Generically, we consider an integral of the type: 

»+oo r+oo r+ca 

dk did duo' 5 {1 — uo — uj') f{k,uo,Lo')pi{k,uj)pj{k,uj'), (4.5) 

J —oo J —CO 

where i, j = t, /. The integral over uj' can be used to eliminate 5 {1 — uj — u'). We get: 

+00 f+oo 

dk / dujf{k,uj,l — uj)pi{k,uj)pj{k,l — uj). (4.6) 



Using ()4.2|) . there are three kinds of contributions, namely 66, 6<d and G0 contribu- 
tions. In all cases for the indices i and j, the 66 contribution is always zero because of 
kinematics. Indeed, in order for it to be nonvanishing, the gluon energies must satisfy 
±uji{k) ± ujj{k) = 1, relations always forbidden by the dispersion relations ()2.H|1 . The two 
6Q contributions are worked out as follows. For example, we want to get an expression 
for Jl^'^dk J^^dujf{k,uj, l-uj)^i{k)/3j (fc, 1~uj)[6{uj~ uji) - 6 {uj + uji)] Q {k - |1 - uj\). A 
non-zero contribution must satisfy u = ±uJi{k), together with 1 — k<uj<l + k. It 
is not difficult to see that only the case u = uJi{k) is allowed, and this for all values 
k > p. The integration over uj is straightforward and we obtain for our term the result 
J^'^dkf{k,uJi, 1 — uji)}if3j {k, 1 — tOi). The other 6Q contribution is obtained in a similar 
manner and is equal to j^"^ dkf{k, 1 — ujj, uj) ij jSi {k, 1 — Uj). The 99 contribution does 
not vanish when u satisfies simultaneously —k < uj < k and 1 — k<uj<l + k. This 
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amounts to having k > 1/2 and 1 — k < u < k. We write then the 00 contribution 
Ii/T^^ Ii-k^^f(^^^^ ^ ~^ 1 ~ Putting all these intermediary results 

together, we have: 

dk [f{k, uji, 1 - uji) 3i pj {k, 1 - cjj) + f{k, 1 - ujj, Uj) pi {k, 1 - Uj)] 

p+oo pk 

+ dk dujf{k,uj,l-uj)pi{k,uj)pj{k,l-uj). (4.7) 

Jl/2 Jl~k 

Another accompanying type of integrals that can straightforwardly be deduced from Iq 
is the following: 

/■+00 z'+oo r+oo 

I^= dk dcu duj'5{l-uj~u;')f{k,ij,uj')e{k-\uj\)pj{k,uj') 

J fj, J —oo J —oo 

"+00 r+oo 

dk I dujf{k^uj,l — uj)Q{k—\uj\)pj{k,uj') 

' fl J —oo 



j'+oo r+oo pk 

/ dkf{k,l — ujj,ujj) }j + / dk 

Jri Jl/2 Jl-k 



dujf{k,uj,l-uj)pj{k,l-uj). (4.8) 

' fl Jl/2 Jl-k 

Indeed, it suffices to replace in ()4.7|) 3j(A;) by zero and Pj{k, 1 — u) hj one. 

Only integrals of type Iq and Jq contribute to the coefficient of zeroth order in p in the 
damping rate. With the appropriate functions /, the remaining integrals in ()4.7j) and ()4.8|1 
are all finite in the infrared, and so fi can safely be taken to zero. However, it is important 
to mention that these two types of integrals Jq and Iq do also intervene in the second-order 
coefficient in p and, with the corresponding functions /, are infrared sensitive. We will see 
this explicitly later in the text. 

The second type of integrals we have to deal with is the following: 

P + OO P + OO P + OO 

h = dk duj j duj' 5 {1 — uj — uj') f{k,uj,Lu')pi{k,uj)dkPj{k,uj') 

J pL J —OO J —OO 

OO P + OO 

dk / dujf{k,uj,l — Lu)pi{k,Lu)dkPj{k,l — uj). (4.9) 

fl J —OO 

Here too the discussion has to be carried out contribution by contribution, using the structure 
of the spectral functions ()4.2j) . The first contribution to consider is the one that involves two 
delta functions, for example j^'^dk j^^duo f 5 {uj — Ui) dk [^jS {1 — u — Uj)]. When the 
derivative over k is applied to ij{k), the supports of the two delta functions do not intersect 
and we have zero contribution. But this is also true when dk is applied to 6 {1 — u — ujj{k)) . 
This is because whatever the result of the action of dk, the two delta functions do not vanish 
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only on the intersection of their respective supports, which is always empty by kinematics. 
This can be checked explicitly by regularizing either the derivative over k (i.e., replacing 
it by a finite difference) or the distribution 5 {1 — cu — Uj{k)). In fact, the same argument 
applies when we have a second-order derivative over k and, in general, terms involving 
6 {u ± Ui) d]^ [6 {1 — CO ± ujj)] do not contribute. 

The second contribution in Ii to consider is an integral involving Sdk (/?©). Consider for 
example J^°°dk J^^dujf{k, uj,l - uj) {uj T ^^i) dk \j3j {k,l - uj) Q{k - \1 - uj\)]. Applying 
dk to Pj is not problematic: only the case cu — uJi{k) with k > /i contributes and we obtain 
for this f^°°dkf{k,uji, 1 - Ui) li dkpj (fc, 1 - Applying dk to e(A; - |1 - uj\) gives 

6{k—\l—uj\) and the support of 6 {u ^ uji) 6 {k — \1 — uj\) is = 1 and k = 0. Since 
k > > 0, this part of the contribution is zero. 

There is another contribution that involves a S and a 6, namely J^°°dk J^^du!f{k, uj,! — 
uj) Pi {k, u;)Q{k — \u;\) dk {^j S {1 — lo^ ^j))- We first apply the derivative over k to and 
obtain the piece j^°°dkf{k^ 1 — {k, 1 — cuj), where 3^- stands for dij{k)/dk. We 

then apply it to 5 {1 — u ^ ujj). We use the standard rules regulating the handling of 
the delta distribution and we always check the results by regularizing either the deriva- 
tive dk or the delta function itself. Only uj = 1 — ujj contributes and we get the piece 
J^'^dku'^lj [f{k,l-u, u) (3i {k, 1 - u;)] |^^^ . . 

There is one last contribution that we have to consider in Ii, namely the double integral 
J^°°dkj;^^du}f{k,u, l-ou)f3i {k,u)e{k - \cu\) dk [f3j {k, l-u)e{k-\l- The appli- 



cation of dk to Pj {k, 1 — cu) gives Ji/^dk J-^_j^dujf{k, a;, 1 — a;) (3i (k, cu) dk(3j (k, 1 — a;) and its 
application to 6 (A; — |1 — cj|) gives zero because of kinematics, as before. Putting all the 
above contributions together, we obtain the result: 
/•+00 

Ii= dk f{k,u}i,l-u;i)iidkPj{k,l-Lj)\^^^,+ f{k,l-LJj,LJj)i'j(3i{k,l-ujj) 

J u 



+^ Ij [f{k, l-u;,u;) (3i (k, 1 



cu 



p+oo pk 

+ dk dcuf{k,cu,l-cu)pi{k,cu)dkpj{k,l-cu). (4.10) 

Jl/2 Jl-k 

Like for /q, there is here too an accompanying integral I[ for Ii which can be deduced 
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straightforwardly from ()4.1Up by replacing 3j by zero and Pi {k^oj) by one: 

;"+oo r+oo r+oo 

I[ = dk duj I duj' 5 {1 — uj — uj') f{k,uj,uj')Q{k — \uj\)dkPj{k,uj') 

J fi J —oo J ~oo 

oo r+oo 

dk / dujf{k,iu,l — uj)Q{k—\uj\)dkPj{k,l—uj) 

fi J — oo 



-oo 



dk 



f{k,l-Uj,ujj)^'j +^3j d^[fik,l-uj,uj)]\^ 



r+oo rk 

+ dk dujf{k,iu,l - uj)dk(3j {k,l - uj) . (4.11) 

Jl/2 Jl-k 

The third type of integrals we have to deal with is one that involves a second derivative 
in k: 

+00 r+oo r+oo 

dk dijj I duj' 6 {1 — uj — uj') f{k,uj,uj')pi{k,uj)dlpj{k,uj') 

fi J —oo J —oo 

-oo r+oo 

dk dujf{k,uj,l-uj)pi{k,uj)dlpj{k,l-uj). (4.12) 

fi J —oo 

The steps to treat the different contributions parallel those followed for Ji. As explained 
before, the 66 contribution is zero because of kinematics. Therefore, the first contribu- 



tion to look at is °^dk J_^du!f{k, to , 1 — to) 6 (to ^ cui) (9| [Pj {k,l — to) Q(k — \1 — tu\)]. 
We know that dl {PjQ) = dlpj 6 + 2dkPj dk& + Pj 9^6. The first term dlpj 6 yields sim- 
ply f^'^dkf^kjtOi,! — tUi) dlPj {k,l — tu)]^^^ . The second term 2dkPj dkO yields zero 
because the kinematics imposes k = , which is excluded from the integration region 
as explained. The same is true for the term Pjd^Q = Pjdk6{k— |1 — c^l) because here 
too the support of the delta functions forces uj = 1 and A; = 0, which is excluded. The 
next contribution is f^°°dk f^^dujf{k, u,l — to) Pi {k, to)Q{k — \to\) dl [3^ 5 (1 — =F ^j)]- 
We have dl{ij6) = dj 6 + i'jdk6 + }jdl6, where stands for d^^j/dk'^. The term 
l'- 6 yields J^°°dkf{k,l — tOj,Uj) ^'jPi{k,l — toj). The two other terms i'jdk6 + 2ij ^l^ 
yield together [{^,to] + ^J^) + ^.uofdl] [f{k, l-to,to) A {k, 1 - to)] l^^^. The 

procedure is to replace dk with to'^d^ and apply the usual rules regulating the han- 
dling of the delta distribution. Also, integrands that force upon us the condition 
k = are systematically excluded. It remains to look at J^°°dk J^^dtof{k,to,l — 
uo) pi {k,to)e{k - \to\) dl [pj {k, l-to)e{k-\l- to\)]. The first term dlPj 6 yields the 
piece Jy'^dk J^_f^dtof{k,to,l — to) Pi{k,to) dlPj {k,l — to) and the second term 2dkPj 6 
the piece 2 Jy'^dkf{k,l — k, k) Pi {k,l — k) dkPj {k,u)\^^f^, both worked out as previ- 
ous similar terms. The third term in this contribution j^'^dk J^^dtof PiPjQdk6 = 
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C dk Cyuj [dk if PiPjQS) - dk if A) PjQS - f PidkPjQS - f PiPjSS] . The support of the 
Q5 distribution is a; = 1 — k and k > 1/2 and that of the 55 distribution the point 
uj = k = 1/2. Because Pj{k, k) = for every k, the terms —dk (/ A) PjQ5 — f /3i/3j55 are 
identically equal to zero. The term —f PidkPjQS gives the contribution — fy'^dkf{k,l — 
k, k) Pi {k, 1 — k) dkPj {k,Lj)\^^j^. Generally, the term dk (/ PiPjQS) would be dismissed as 
a total derivative. This would be straightforward if we are allowed to integrate first over 
k than over cu. But the permutation of the two integrations is not always permissible and 
care must be taken; each case has to be treated individually. Here, one can show that this 
term gives indeed zero contribution by either regularizing the delta function or replacing the 
derivative by a finite difference. For example, if we use the second method, we would have: 

f+oo f+oo 

/ dk du;dk[g{k,uj)e{k-\uj\)5{k-\l-uj\)] 

J II J —oo 

f+oo f+oo 

c^l/e / dk / du;[g{k',u;)e{k'-\u;\)5{k'-\l-u;\)],,^,^^ 

J fj, J —oo 

f+oo f+oo 

-l/e I dk 1 dujg{k,uj)Q{k-\uj\)5{k-\l-uj\) 

J n J —oo 

where e — >^ and g{k,uj) stands for f (5i(5j. This difference is equal to l/e fy'^dk'g{k' , 1 — 
k') — l/e J^,'^dkg{k, 1 — k) = 0. Putting all the above results together, we obtain: 



'1/2 
r+oo 
/ dk 



f{k, uji, 1 - Ui) li dlpj {k, 1-Lj) 1^^^. + f{k, 1 - ujj, ujj) I'-pi {k, 1 - ujj) 
+ [{iju'; + 23;^.) d^ + ijufdl] [f(k, l-u;,u;) A (k, 1 - u)] |^_." 



r+oo 

+ / dkf{k,l - k,k) Pi{k,l - k) dkPj {k,^)L=k 
J 1/2 

r+oo fk 

+ dk dujf{k,uj,l - u) pi{k,uj) dl(3j {k,l - u) . (4.13) 

Jl/2 Jl-k 

As for the two types Iq and Ii, we write for the accompanying integral 

f+oo f+oo f+oo 

72=/ dk dcu duj'5{l-uj-cu')f{k,uj,uj')e{k-\cu\)dlpj{k,(jj') 
J Ii J —00 J —00 

f+oo f+oo 

— dk dLjf{k,LJ,l — Lj)Q{k—\u}\)dlpj{k,l — Lj). 
J II J —00 

r+oo 

= / dk f{k,i-ujj,ujj)3';+ {{i,u; + 2^yj)d^ + ^,ufdi)f{k,i-uj,u)i^^ 
J II 

f+oo r+oo fk 

+ / dkf{k,l-k,k)dkl3j{k,uj)\^^j^+ j dk j duf{k,uj,l-uj)dli3j{k,l-u) 

Jl/2 Jl/2 Jl-k 



(4.14) 
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Before carrying through, it is important to note that whereas f3j{k,k) = for all k, the 
partial derivative dkl3j{k,uj)\^^i^ present in ()4.13l ) and ()4.14|) is infinite. This divergence 
comes from light-cone momenta and a more appropriate notation would be dk(3j{k,uj)\^^f^. 
Such terms would spoil the calculation if we were not fortunate to have them cancel with 
similar terms hidden in the last contribution in (|4.13p and (j4.14p . This will be shown later. 
The next type of integrals we must handle is: 

dk duj I duj' {1 — uj — uj') f{k,uj,uj')pi{k,uj)pj{k,uj') 



J —OO J —OO 

D n+oo 

dk I du d^[f\k,uj,J)pi{k,uj)]\^,^-^^^pj{k,l-uj). (4.15) 

J — OO 

Here too the residue-residue contribution is zero. The residue-cut contribution is equal to 
- j^'^dk j^^duiid^ [f{k, u, uj')S (u =F uji) Pj (k, u') Q{k- \uj'\)] \^,^-^_^. Actually, the easi- 
est way to handle this contribution is to reconsider in ()4.15|) a derivative with respect to 
u' instead of u, i.e., replace d^6 (1 — u; — u') by {1 — uj — oj'). Then, the residue-cut 
contribution would write: 



hoo 



dk I duii5 {uj ^ Ui) duji [f{k,uj,uj')(3j {k,uj')Q{k — \uj'\ 



'ijj' = l—U) 



+ 00 



dkii [f{k, Ui, uj)pj{k, u) 



'UJ = l—UJi 



where we have implicitly used the usual supports of the delta and theta functions. 
The cut-residue contribution is worked out in a similar way and found to be equal 
to — J^°°dk^j duj[f{k,uj,ujj)Pi{k,uj)]\^^-^_^ . The last contribution, the cut-cut term, 
writes - fl^'^dk f^^dud^ [f{k, u, {k, u) {k, u') Q (k - \uj\) Q {k - \uj'\)] |^,^,_^. Ap- 
plying the derivative over uj on 00 implies that either uj = k oi uj' = k. In both cases 
j3(3 = 0, which means no contribution from here. Hence the cut-cut contribution is simply 
equal to - J^/^dk J^_j^duj(3i {k, uj) d^, [f{k, uj, uj')(3j {k, uj')] So we get: 



dk 



OO 



3i [f{k,uji,uj)l3j{k,uj)]\^^^^^^ +}j [f{k,uj,ujj)(3i{k,uj) 



dk / duj(3dk,uj)d^, [f{k,uj,uj')(3,{k,uj')]l,^^_^. (4.16) 

1/2 Jl-k 

The accompanying integral Jg where pi{k,uj) is replaced by Q{k — \uj\) is obtained from ( 
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I4.16P as usual. We have: 



/+00 r+oo 
duj / duj'dJ{l-uj-uj')f{k,uj,uj')Q{k-\uj\)pj{k,uj') 
^_ -oo J —oo 

oo r+oo 

dk I duj d^[f{k,uj,uj')Q{k-\uj\)]\^,^^_^pj{k,l-uj) 

fi J — oo 

oo r+oo rk 

dk^j d^f{k,uj,ujj)\^^ - dk / dujd^, [f{k,u;,uj')Pj {k,uj')]\^,^^_^ . 

H Jl/2 Jl-k 

(4.17) 

The next type of integrals we have to deal with is: 



h= I dk du duj'5{l-LU-Uj')f{k,UJ,Uj')6{k^-UJ^)pj{k,Uj') 

J —oo J —oo 

x> r+oo 

dk / dujf{k,uj,l-uj)S{k'^ -uj^)pj{k,l-uj). (4.18) 



' /i J —oo 

It is not difficult to work it out. Using ()4.2|) . the delta-delta contribution has zero support 
and hence vanishes. The delta-theta contribution is easily worked out along the usual lines 
and we obtain: 

dk 

^4= / ^f{k,k,l-k)P,{k,l-k). (4.19) 

There remains one last type of integrals to look at and which too is not difficult to 
manipulate. It is: 

r+oo r+oo r+oo 

h= dk duj duj' 5 {1 — UJ — u') f{k,uj,uj')duj[6{k — uj) — 5{k + uj)] pj{k,uj') 



fj. -J — oo -J — oo 

oo r+oo 



dk dujf{k,uj,l-uj)d^[6{k-uj) -6{k + uj)]pj{k,l-uj). (4.20) 

' fi J — oo 

Here too the delta-delta contribution is zero because of kinematics. The delta-theta contri- 
bution is treated by turning the derivative onto fjSjQ. Manipulations by now familiar yield: 

30 

dk [fik, LU, 1 - uj)^,ik, 1 - cu)]|_, . (4.21) 

'1/2 

Every integral in the expressions (IT^ . (jSHSl), (E2ZI), and (ESI) falls into one 

of these above types. We are therefore ready to perform them. However, as we mentioned 
more than once in the text, there is occurrence of infrared sensitivity to be aware of and 
which we need to treat with care. There are actually two types of divergences: infrared and 
light-cone. The ffist type will persist in the final result and will be expressed in terms of the 
magnetic mass p we introduced from the start as a physical infrared regulator. The second 
type will not persist as we will see in the sequel. 
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B. Infrared and light-cone behavior 



It is best to work out few sample examples in some detail, which all come from the 
coefficient of in the sum (j4.ip . Let us start with the pip'i contribution. Up to the factor 
g'^NcT/ 12071 for which there is no need writing explicitly, the corresponding / function is 
f{k,uj,uj') = 3A;2 (3 + AAk"^ - 12ujuj') /2ujuj'. This integral is of type Jq. From (gTj), we can 
expect infrared divergences coming from the two terms, which in fact are equal because 
here f{k,uj,uj') is symmetric in u and u'. Using the small-/c expansions of on-shell energies, 
the residue and cut functions given in appendix A, we have the following small-/c behavior: 



f {k,l-uji,uji)^i(3i{k,l-uji) = ^ + 

Sk 



2217 



27 



400 
10230949 



-TT 



3200 
134397 



k 



243 



-TT 



-TT 



k^ + 0{k^), (4.22) 



1680000 1120000 1280000 

which indicates a logarithmic divergence. The corresponding integral can be written as 
follows: 

r> + 00 



6 dk 6 , , , 
8Xt + 8'"*-W 



2 / dkf {k, 1 - uji, uji) ^i{k)(3i{k, 1 - ui) 

+ 2 dk 
Jo 

•1 

+ 2 1 dxk[{x) f{k,l-ujuuji)ii{k)(3i{k,l-uji\^^^^^^ (4.23) 



f {k,l- ui,ui)2ii{k)Pi{k, 1-uJi) - 



8k 



What we did is this. We split the original integral in two: one from /i to a finite value plus 
one other from that finite value to +oo. The second integral is finite, and it is (numerically) 
better for it to change the integration variable from A; to x = k/uji{k), which implies that 
k = ki{x). The finite value in question that splits the original integral is chosen in terms of 
X instead of k and is denoted by i. In the integral from fi to ki{i), we subtract Q/8k from the 
integrand, which renders it safe in the infrared and so we can take the lower boundary /i — ^ 0. 
We must add of course the contribution | /^^'^^^ f = | /^^ f + I ^^(^)- ^he finite value i 
must be chosen small enough in order to make the integral J^''^^^ dk [f^iPi — ^) numerically 
feasible. Indeed, though we are assured of its finiteness analytically, both integrands f^iPi 
and ^ are still each (here logarithmically) divergent. However, for a small value of i, say 
i = 0.1, we can use an expansion of [f^iPi — ^) in powers of k in order to get a number for 
the integral. We have pushed the expansion to 0{k^^) and we start having good convergence 
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for already £ = 0.6. Also, we must (and do) check that the final result does not depend on 
a particular value of i. We finally get: 

2 dkf {k,l-uji,uji)^i{k)(3i{k,l-uji) = -^\n 12 + 5A6353171. (4.24) 

It is important to note that the infrared divergence in ()4.24p is not coming from an 
eventual expansion around ultrasoft p of the dressed propagator *A(P — K) or any similar 
distribution. Actually, as mentioned in the paragraph after ()4.8|) . a pip'i term is present in the 



zeroth-order coefficient of the damping rate 



48l ]. but with a corresponding function / going 



sufficiently fast to zero so that the divergent behavior of ii{k)(3i{k, l—ui,) ~ —l/AQk + 0{k) 
is neutralized. What happens here is that our function / (fc, 1 — ojuuji) ~ —15 + 0(/c^), a 
behavior not fast enough to screen the infrared divergence. Therefore, it is clear that it 
is not the expansion of quantities similar to the dressed propagators *A(P — K) around 
ultrasoft p that is sole responsible for the occurrence of infrared divergences as one might 
think. Note also that the same thing happens to the other contributions to the imaginary 
part of the one-loop HTL-dressed self-energy (t/, tt, etc.). 

Returning to the integral we are handling, the corresponding double integral is 
j^l^ dk dujf{k, uj,l — uj)(3i{k, uj)(3i{k, l~uj) = 2.458381509 with no problem worth men- 
tioning. We finally get: 

V- (3 + 44A;2 - Uujuj') pip'i = --\np + 7.921913219 , (4.25) 



where f T> stands for f °°dk f °°dtu f °°dco'5 (1 — cu — u') for short. 

Along the same lines and using the type Jq given in ()4.8|) . we get the accompanying result: 



f) - 4) = "i^^^ + 0.611143948. (4.26) 
In terms of the foregoing discussion, convergence is already obtained for I = 0.5. Note that 
here too, there is a (logarithmic) infrared divergence without involving an expansion of the 
dressed propagators and similar distributions. 

The next generic contribution in ()4.Hl we consider is the one that involves pidkp'i, with 
a corresponding function f{k,uj,u!') = Ak^ {13 + Sk"^ + lOu — du"^) juJio'. The type of in- 
tegral we need to use is l\ given in ()4.1()j) . Let us look at the term f^idkPi- Using the 
expression of / and the small-fc expansions of and d^Pi given in appendix A, we have 
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fM - -i - (1 - If ) + - ^ - ii;) + Oik% hence a logarithmic 

divergence. We operate along similar lines as indicated in the case pip'i above and ob- 
tain dkii{k)f{k,uji,l - uji) dkPiil = f In^ - 0.404833413. The conver- 
gence is good starting from £ = 0.3. Next we look at the term i'if(3i. Its expansion 
in powers of small k yields I'JPi — ^ -\- + ^f ^ k -\- 0{k^). A logarithmic diver- 
gence too. We get dk^'ifPi — yln/^-l- 2.193663437, with a good convergence starting 
from i — 0.5. Next is to look at co'i^idi^ {f Pi) . We have the expansion u'liid^ {f Pi) — 
- (i + If ) {§M + ^iif + iiS) + This integrand converges fast enough 
so that no infrared singularity is present and we have f^^dkuj'i^id^ if Pi) = —2.222237940. 
This finite result also means that not all contributions are divergent in the infrared, partic- 
ularly regarding those coming from an expansion of the effective propagators and indeed, 
we encounter many such instances. 

We look now at the double integral in pidkp'i, namely dk J^_^ doi)A(k, u) where A(k, u) 
stands for PifdkPi- Here care must be taken because the integrand, more precisely dkPi, is 
divergent when u 1 — k. To handle this, let us make the change of variable y = uj — 1 + k 
in place of lu and the change of function G{k, y,Y)/y = A{k, 1 — k + y), where Y = 1/ \ny, 
small when y is small. It is easy to check that G{k,y,Y) has no pole in y. Make then the 

expansion G{k, y, Y) = Go(^, Y) + Gi{k, Y)y + O (y^). We obtain the result: 

27A; (56 - 24^;^ + 32A;) 



Go{k,Y) 



4 ((3 + + ^ In (2A; - 1))' + f^^^) 



(|(Fln2A;-l)-(A;2 + 3)F)F3 

X ^ ^ ^ ^—^ 2- (4-27) 

(((3 + F)y - I (y ln2A; - 1))' + ^n^Y^^ 

Gi{k, Y) has a lot more comphcated expression that we do not need to display. A small Y 
expansion of GQ{k,Y) yields: 

= -2(56.3 24^. ,.,., ^- + 0(n. (4.28) 

Therefore, GQ{k,Y) is smooth and finite when F ^ and the divergence of A{k,uj k) 
goes only like 1/y. This suggests to rewrite: 

/•oo i-k /•oo /■2fe-l 7 

dk dujA{k,uj)= dk —[G{k,y,ln-^y)-Go{k,ln-^y)] 

Jl/2 Jl-k Jo.5+ Jo+ y 

/•oo /•o- 

+ dk Go{k,Y). (4.29) 

^0.5+ 7l/ln(2ifc-l) ^ 
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Subtracting Go{k,l/ \ny) from G {k,y,l/ \ny) in the first integral makes the integrand 
finite and y can safely go to zero. It is Go/y itself that is integrated in the second term, 
but the integration variable is changed from ytoY = l/\ny. Since G'o(^)^) starts as Y^, 
see (j4.28p . the limit Y —>■ becomes safe. Therefore, the light-cone divergence present in 
A{k,uj) = j3ifdkl3i is circumvented by the above adequate change of variable y ^Y. The 
first integral gives —2.0291277 with a very good convergence at fcmax = 100 and 0+ = 10~^^. 
The second integral is equal to 0.253206577. Therefore: 

oo nk 

dk rfwA (A:, /(fc, 1 - cj)9fcA (fc, 1-cj) = -1.775921122. (4.30) 

1/2 Jl-k 

Putting all the above partial results together, we arrive at: 

/ V (13 + 3k^ + lOu - 9u;2) pi dtp', = — In - 2.209329038 . (4.31) 

J ujuj' ^ ' 3 

Before tackling the next example, two comments are in order. First, as we see, light-cone 
divergences are brought under control by a suitable change of variable whereas the infrared 
divergences are not. This situation is the same with all the remaining terms. Furthermore, 
as we see, the handling of the light-cone divergences does not necessitate the introduction 
of an asymptotic gluon mass moo as in the improved HTL summation proposed in jsS^ . 
However, we should mention that we do not know if this works for other physical quantities 
as we do not have a general proof for it. 

The second comment is about A^max = 100 we mentioned as the upper bound for k for 
which the (numerical) integration starts to converge satisfactorily. Remember that all is in 
units of m^, the soft scale. Remember also that at dressed one loop, all energies and momenta 
have to be soft. One asks then: is lOOm^ still soft? This point is pertinent because most 
similar integrations necessitate a fcmax in the hundreds. This observation is an indication 
that for the HTL approximation to be applicable, we have to be well in the quark- gluon 
plasma phase, i.e. at quite high temperatures so that lOOm^ can still be considered soft. 



1 _L ±i£ do; 



The 



The next contribution we discuss is 6(9fcP with f{k,u!,uj') = ^ 
corresponding integral is I[ given in ()4.11|) . The first contribution we consider is /j^. A 
small-fc expansion yields /j^ = 9/k^ + — ^^k + Oijc"). Here we should expect a 1/ p? 
behavior plus the usual In p. We therefore write dkfi'i = ^ — 2k^j)2 — tm 



^\nki{£) + Jq dk [f^'i - ^ - j^] + dx},'J\k=ki(x)- Now the integrations are smooth 
and good convergence and £— independence start at £ = 0.6. We obtain J^°° dkfi'i = 
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— In /i — 1.051103551. The next contribution to look at is uj'i^id^f. A small-fc expansion 
yields uj[^id^f = ^^ + ^k+0{k^), only logarithmic a divergence. Using the same technique, 
we get dkuj'i^idujf = || In /i + 0.3027125835 with good convergence. 

Here too the corresponding double integral has a divergent integrand as 1 — k, 

but the light-cone divergence can be circumvented by a simple integration by part. Write 
d{k, u) = fik, l-u,u) = [3 - ^ + (15 - ^) - (12 -^)uj^- ^uj^] . We have: 

oo pk poo pk 

dk I dujd{k,uj)dkPi {k,uj) = — / dk dujdkd{k,uj) (3i{k,uj) 

1/2 Jl-k Jl/2 Jl-k 

oo pk ^ 

dkd{k^l — k) I3i{k^l — k) +Yim I dtod{k,uj) Pi{k,uj) 

1/2 Jl-k k^l/2 



(4.32) 



The first integral is finite. Indeed, write it as Jq dK J^^/J^^ ^^^^^^(I/k, cu) /?/ (I/k, cu) where 
K = 1/k and it yields —0.2338436626. The second one is also finite and rewritten as 
— ^d{l/ K,l — 1/ k) Pi (1/ K,l — 1/ k), it yields exactly —0.24. Note that we have no 
independent analytic check of this exact value. The third term can be rewritten as 

1/k k^O 

lim dudil/ K,u!) Pi (1/ K,u) and the integral can easily be shown to go to zero 

''^ K^2 

in both limits n and n ^ 2. The light-cone divergence is thus circumvented. Putting 
all the above partial results together, we get: 



V — 



9 27 

eafcp; = — + In u- 1.22223463. (4.33) 

'^^ 2/i2 100 ^ ' 



The next contribution we consider is pidlpi. The corresponding integral is of type 
I2 given in ()4.13j) and the function f{k,LU,uj') = ^(2 +5a;). The treatment is as be- 
fore. The first term to consider is I'lfPi- A small-fc expansion yields I'l f Pi = f + 
^— ||i — k + 0{k^)^ hence a logarithmic divergence. Using the same techniques as 
above, we get dk^i'lfPi = — 21n/i — 1.435815399, with a very good stability starting 
from i = 0.5. The next term to consider is dk^ifdlPi. The small-A; expansion yields 
hfdlPi = 4 + (i - w) ^ + 0{e). We obtain /~ dk^JdlPi = -| In/i + 0.731445450, 
here too with very good stability starting from i = 0.5. Next we look at 2 dku'i^'id^{fPi). 
The small-/c expansion yields 2u'ii'idi^{f Pi) = + k + 0{k^). This term is finite and 
we have 2 dkuj'ii'idi^[f Pi) = 1.347361846, with a good convergence. The next contribu- 
tion is 3/ (c^pS^ + ^^I'du,) if Pi)- A small-fc expansion yields — k + 0{k^), finite. We 
therefore have dk^i {jjj'j^d'^ + oj'ld^^) (fPi) = —0.1176374160, with good convergence too. 
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It remains to handle Jy^^^f ^i^kPi + J^2^^ Ji^k^^f f^^'^kPi- We have aheady remarked 
that dkPi {k,uj)\^^i^ is infinite. First a suitable change of variable from u to 1 — u. Then 
perform the second integral by parts to get: 

POO pk pk fc— >oo 

dkfPidk/3i+ dk dufl3idl(3i= dudik,tu)dk(3iik,u) 

1/2 Jl/2 Jl-k Jl-k k^0,5 

POO POO pk 

- / dkdik, 1 - k) dkl3i (k, u) l^^^^ - / dk diudkd{k, uj)dkpi (k, u) , (4.34) 

Jl/2 Jl/2 Jl-k 

where we used the definition d{k, u) = f{k, 1 — u, uj)l3i {k, 1 — u). The first term is equal to 
zero (the integral over as a function of k vanishes in the two limits k oo and k 1/2). 
The second term is also zero because jSi {k, k) = 0. The integrand of the third term diverges 
both in the limit uo k and c<j — > 1 — /c. This can be handled in the following manner. 
Define A{k,u) = dkd{k,uj)dkl3i {k,uj); Fi{k, y) = yA{k,k-y) and F2{k, y) = yA{k,l-k + y). 
Define for i = 1,2 the function Gi{k,y,Y) = Fi{k,y) with Y = 1/lny. An expansion in 
powers of small y yields Gi{k,y,Y) = Gio{k,Y) + O (y) and an expansion of Gio{k,Y) in 
small Y yields Gio{k,Y) = Ki{k)Y^ + 0{Y^) and G2o{k,Y) = K2{k)Y^ + 0{Y^) with Ki{k) 
a (complicated) function of k that does not need to be displayed. Therefore, we can write: 

POO pk POO p(k-0.5) J POO p(k-0.5) j 

- dk dujA{k,uj) = - dk ^Fi{k,y)- dk JiF2{k,y), 

Jl/2 Jl-k Jo.5+ Jo+ y Jo.5+ Jo+ y 

(4.35) 

with: 

POO pk-O.b J POO pk-0.5 J 

- dk JiF,{k,y) = - dk Ji[F^(k,y)-G,o{k,ln-'y)] 

Jo.b+ Jo+ y J0.5+ Jo+ y 

°° /"°~ dY 

dk / G,o{k,Y). (4.36) 

0.5+ ^l/ln{fc~0.5) ^ 

Both integrals {i = 1,2) are finite and the logic behind this technique has been discussed 
previously, around dOHl)- We obtain -0.0059819227 for i = 1 and 0.917793484 for i = 2. 
Hence, in all we get jy^dkf(3idk(3i + J^^dk /^_^ rfw/A^^A = 0.911811562. Putting all these 
partial results together, we get: 

/2k^ 13 
V (2 +5uj) pidlp'i = In /i + 1.437166043. (4.37) 
uu' 3 

The accompanying term to work out is Qdlpi. The corresponding function is / {k, uj, uj') = 
and the type of integral we need is given in ()4.14j) . The first term to look at is /jj'. 
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A small-fc expansion yields fi'l = + — ^^k + 0{k^), hence a plus a logarithmic 
divergence. Using the same techniques, we obtain J^°° dkfi'l = — ^ — ^ In /i + 3.32815131. 
The other contributions are worked out similarly and yield only a logarithmic divergence. 
We have J^^^^dk {{^ilu'/ + 2}[iul) + iiuj'i^dl) / = ^ In/i + 0.1270638712. Next is to look at 
the double integral. We have: 

fk POO rk k^oo 



dkf{k,l-k,k)dkPl 
fc^O.5 J 1/2 



pK poo pK, 

dk / duofdll3i + / dkfdkl3i = / dufdkl3i 

Jl-k Jl/2 Jl~k 

poo pk 

- dk diudJdkPi. (4.38) 

Jl/2 Jl-k 

The first term is zero. The second term is analytically finite but must numerically be handled 
with care because it does not converge easily in the limit k 0.5. A change of variable 
k l//c, usually useful, is not of any help here. The best way we find to treat this is 
to introduce a small e such that the integral /q^^^ dk carries with no problem, and handle 
Jq '^^"^ dk using an expansion of the integrand. The first term in the expansion has been 
sufficient. Thus we find this second term (exactly) equal to —0.645. Note that here too, 
this simple and exact number is obtained numerically and has not been shown analytically. 
The third term has a light-cone divergence and is handled in the {y, Y) variables. We find 

it equal to —0.936156340. Putting everything together, we get: 

/3lj2 q 27 

^ W®^^'^^ = ']r^' — In /i+ 1.874058841. (4.39) 



The next term we work out explicitly is pip'id^S with a function f {k,u!,uj') - 
The corresponding integral is of the type J3, given in ()4.15p . We consider first the con- 
tribution —23^(9^ if Pi), the factor of 2 coming from the fact that two terms in ()4.15p 
are equal. The small- /c expansion yields -y + ^ + 0{k^), which means no in- 

frared divergent behavior analytically, but practically, the integration does not converge 
easily when k —>■ 0. We then split the integral into two pieces, one piece equal to 
J^^'W dk [(^-M + j A; + (iff - »^ - ^) k' + Oik')] where we perform an ex- 
pansion of the integrand, and a second piece equal to —2 dxk'i^id^ if Pi) which is integrated 
smoothly. We get in all 0.00866285165, with very good stability and £— independence al- 
ready from i = 0.2. The corresponding double integral presents a light-cone divergence that 
is handled in the {y,Y) variables. We find it equal to —0.825017233. In all, we have: 

V pip'idj (l-u-uj') = -0.8163543814. (4.40) 
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The last term we discuss is Qp'id^S with a vanishing corresponding /. Hence: 

Vf {k, to, uj') Qp[dJ {l-uj-uj') = 0. (4.41) 



These are all the sample terms contributing to the second coefficient in of Im*!!^ we 
wanted to discuss in the main text. All the remaining contributions are treated along parallel 
lines and no new subtleties not touched upon here do occur. The only new feature is the 
occurrence of divergences of order l//i^. We give in appendix B all the partial results for 
each contribution. 

The coefficient of order zero in Im*!!^ is finite and is already determined in Putting 
everything together and restoring back m^, we obtain: 



Im*nt(P) 



^13.086 + ^90.064//i^ + 2.7 /f + ^ In + 1.5192728^ + O (p^) 

(4.42) 

where Ji = fi/rrig, a dimensionless parameter, and p = p/mg. Using (|2.7|) , we finally find: 

7t(p) = ^ [6.543 + (45.032/^^ - 314.21/ fi^ - 507.96 In/i + 1014.65) p"^ + O (p^)] . 

(4.43) 

As we have claimed, when the damping rate 7t(p) for transverse gluons is calculated in 
the sole context of HTL-summed perturbation theory in an expansion in powers of p^, the 
second coefficient is divergent in the infrared. We discuss this result in the following and 
last section. 



V. DISCUSSION 



Using the hard-thermal-loop perturbative framework, we have determined the damping 
rate 7t(p) for ultrasoft transverse gluons in high-temperature QCD to its lowest order g'^T. 
The calculation requires the obtainment of an expression for the imaginary part of the fully 
HTL-dressed one-loop-order gluon self-energy. We have determined this quantity in the 
form of a series in powers of the external momentum p = p/mg, a small parameter when p 
is ultrasoft. The aim of this section is to carry a critical discussion of the results (14.421 ) and 
fl4.43|) . and the method used to obtain them. 
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A. Magnetic mass as an infrared regulator 



The presence of and In/i (to a lesser degree) in the coefficient of in ()4.43l ) 

is troublesome to some extent because it suggests a divergence when /i 0. However, we 
may regard the issue from a different perspective. Let us suppose that the HTL-summed 
perturbation is all of massless QCD at high temperature. Since, on physical grounds, we 
must have 7^(0) = 7t(0), we can impose this condition on { II. Sj) and obtain thus a (crude) 
estimate of the magnetic mass: 

/X = 0.32259 mg + ... . (5.1) 
We can put this value back in (j4.43p and obtain the following expression for 7^ = 7^/777,^: 

7t(p) = [6.543 + 2728.3 p'' + O [f)] g + ..., (5.2) 
where q = f^" q. Note that q is a better expansion parameter than q because it is 

^ 2A-Ky/Nc+Nf/2^ ^ 1-1- ^ 

much smaller. Indeed, for Nc = "i and Nf = 1, we have g = 0.0638 (yf, and for Nc = 3 and 
Nf = 3,g = 0.0563 g. 

Result ()5.2j) is quite acceptable, even with a relatively large coefficient for p^. Remember 
that in massless high-temperature QCD, there is a hierarchy of scales g"'T, with n a nonneg- 
ative integer j^. Consider for instance the quantity fi^ (p) = ^It (p) /nig, where Vtt (p) is the 
full (complex) on-shell energy of the ultrasoft transverse gluon. In an expansion in powers 
of the coupling g, we anticipate the following: 

0, ip) = cutip) + [cui'\p) - tjtip)) + (cuf )(p) - Z7P^(P)) + • • • , (5.3) 

where outip) = ^t{p)/^g = l + 3/5p^ — 9/35p^ + ... from ()2.3j) and is of order one (lowest- 
order contribution); ujl.^\p) and 7t(p) are functions of p times g (first-order corrections); 
^t^\p) s-iid "yt'^\p) are functions of p times g^ (second-order corrections), etc. Remember 
that the very definition ()2.6|) of '~ft{p) is based on the feasibility of the expansion ()5.3|) . The 
full complex energy Qt is a solution to the full dispersion relation: 

-n^t+p^-Ut{-tnt,p) = o, (5.4) 

where Ilf is the full transverse-gluon self-energy, a fortiori expected also to admit an expan- 
sion in powers of g: 

Ut = + ) + nf + . . . . (5.5) 
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The logic behind the hard-thermal-loop perturbation in g or equivalently in g can be under- 
stood as follows. In order to get the lowest-order term ut of Qt, one uses only the lowest-order 
term uf'^ = 6Ilt of in the dispersion relation (|5.4j) . The contribution 6Ilt is obtained by 
summing naive-perturbation one-loop graphs with hard internal momenta, more precisely 
momenta ranging from the soft scale gT, or preferably nig, to +00. In order to get the 
first-order corrections cjf^'' and 7t to Clt, one assumes that all of the next-to- leading-order 
contribution is *Ilt and adds it to 6Ilt in fl5.4|) . It is evaluated as a sum of HTL-dressed 
one-loop graphs with soft momenta ranging from g'^T to gT or, more preferably from a 

physical standpoint, from the magnetic scale /i to the electric scale mg. In order to get 

— (2) — (2) ^ (2) 
the second-order corrections and 7^ to fif, one adds 11^ to the dispersion relation, 

assuming that all of it is a set of HTL-dressed two-loop graphs with ultrasoft loop-momenta 

ranging from the scale g^T to the magnetic scale fi. And so on. At the present time, 

there is no known physical process that involves in a natural way the scale g^T and explicit 

calculations in HTL-dressed perturbation can be carried with much hardship only to first 

order. 

Hence, it is quite admissible to use the magnetic mass fi as an infrared regulator in 
HTL-dressed one-loop calculations, and this we did in the present work. The problem that 



arises is how to estimate fj, itsel 



is a non-perturbative effect 
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. The common belief is that chromomagnetic screening 
The somewhat heuristic argument leading to ()5.1|) 



may be useful, but is probably a temporary way out. It assumes that all of massless QCD 
at high temperature is in the HTL-summed perturbation. This is certainly true to lowest 
order, and some would argue that it is also true to next-to-leading order. But there is no 
certainty beyond. Recall that there is no generating functional that yields the whole series 
of HTL perturbation. The only available generating functionals are those generating the 
hard thermal loops only. 

Still, the result ()4.43|) depends strongly on the infrared regulator and for the term in 
to be regarded as a correction to the lowest-order coefficient, we must have p < 0.049 m^. A 
different perspective is to look at the full complex energy ^It (p) in an expansion in powers 
of p = p/rrig. This must be allowed for ultrasoft p on physical grounds. Write then: 

iP) = (g) + ^ (g) p' + {g)f + ... (5.6) 
Where the different vlf' (^), z = 0, 1, . . . can be written in powers of g. For example, [g) 
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is the full zero-momentum complex energy for the transverse gluon (in units of rUg), and we 
have: 

(^) = 1 + [^l'\0)/9 - 6.543*] g + O [f) . (5.7) 

The quantity ujI:^\o) is the next-to-leading order contribution to the gluon thermal mass, and 
assuming that both the longitudinal and transverse gluons are the same at zero momentum, 
we do have a number for it from |^ in the case of a pure gluonic theory (A'^^ = 0). In our 
notation, ul.^\o)/g = — 22.6\/iVc. We also have, with the estimate ()5.1|1 for /x: 

(^) = - + - 2728.3 i\+0 if) , (5.8) 

where ujH^ is the coefficient of in the expansion of ujl:^\p), and 3/5 comes from the 
expansion of u!t{p) = 1 + 3/5p^ — 9/35^*^ + . . . . We do not have an expression for tufj'', 
but for the perturbation to be reliable, we should at least have 2728.3 (7 < 3/5, which means 
g < 0.0002, or ^ < 0.003 for A^'c = 3 and Nf = 1, and g < 0.004 for A^^^ = 3 and Nf = 3. This 
implies that the temperature T must be quite high so that the (running) coupling constant 
be so small (asymptotic freedom). At very high temperature, classical behavior dominates 
and quantum effects are small. It seems therefore that the HTL-summed perturbation may 
be reliable well into the quark-gluon plasma phase, but if we try to lower the temperature, 
this perturbation may become less useful. 

Another interesting rough inference from this 'orders-of-magnitude' discussion is the re- 
lation of the magnetic mass fi to the magnetic scale g^T. Indeed, let us write fi = ag'^T = 
^agmg = 0.32259 nig. This implies that ag = 0.32259^ = 0.012835, and with g < 0.0002, 
we have the lower limit a = n/g'^T > 64.18. This is between one to two orders of magnitude 
large, than the magnetic .cale and add= wate. to the —on that not only staple 
integer powers of the coupling g play a role in this perturbation [1|]. 



B. A Slavnov- Taylor identity 

Away from these orders-of-magnitude estimates, we focus attention on a Slavnov- Taylor 
identity for the gluon polarization tensor derived in in the Coulomb gauge and apphed to 
the next-to-leading-order contribution. The identity recovers already known transversality 
conditions, and when applied to the imaginary part of next-to-leading order gluon self- 
energy, it implies the equality between the damping rates for non-moving transverse and 
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longitudinal gluons. This result is important and worth discussing in the present context. 
We first recall the main steps of that work with a slight change in the notation suitable for 
our purposes. The Slavnov- Taylor identity derived in |55|| is: 

P-n,. (P) = - KP2 - P,P^ + li^ (P)] (P) , (5.9) 

where H^^ (P) is the full gluon self-energy and Ilg;^ (P) is related to the Coulomb-ghost 
self-energy Ilg via the relation Ilg (P) = p^IIgi (P). Note that identity ()5.9|1 is exact and 
valid for all momenta P. 

In the sequel, concentrate on the soft momentum region and consider the expansion of 
the full li^u in powers of the coupling, relation { 15. 5|) . Again, li^^y is the full leading-order 
contribution of order (^^T^, n[jl] the full next-to- leading-order contribution of order g {g'^T'^), 
and so on. When reading the identity ( 15. 9|) in perturbation in powers of the coupling, we 
can deduce the two transversality conditions: 

p^n^) (P) = 0; P'^P'^nW (P) = 0, (5.10) 

both relations known to be satisfied by the hard thermal loop 611 and the HTL-summed 



one-loop-order self-energy *n respectively. For our 
self-energy only to lowest order, and it is derived in 



jurposes, we need the Coulomb-ghost 
to read: 



^'^liP) = -'-iFV^ (5-11) 



which is already of order g'^T and has the limit Ilg^^ (cu, p = 0) = 0. From this we get the 
relation j55| : 

P'^nS (P) = *Ar^ (P) (l, "^p) , (5.12) 

16 p \ P J ^ 

where *Ai is the longitudinal gluon propagator to lowest order given in ()2.1|1 . The above 

relation was checked in to be satisfied by *n as defined in ( 12. 8j) . Since the right-hand 

side of ()5.12p is real, we immediately obtain: 

P'^ImnW (P) = 0. (5.13) 

From here, it is straightforward to show that [ssj ]: 

H{0) = lt{0), (5.14) 
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with the same definition of the damping rates as the one we adopt. One thing pertinent in 
the result ()5.14|) is that, since *n as defined in ()2.8|) is shown exphcitly to satisfy ()5.12|) . it 
too would lead to the same equality, independently of the issue we raise whether *n is the 

full n(i). 

We have emphasized at several occasions that at zero momentum, the longitudinal and 
transverse gluon damping rates have to be equal, and that any consistent computational 
scheme ought to yield this equality. The result of 155| is an additional confirmation of this. 

n nn 

However, all calculations in j48| leading to ()1.H) and in j53|, |5J| leading to (II. Sp . and even 
those of 1^ leading to ()1.2|1 for the quarks are explicit, whereas those in are formal: 
the loop-integral/Matsubara-sum are never done explicitly; no infrared cutoff is introduced 
to regulate the infrared sector and there is no mention of what region should be kept in the 
integral/sum. All these issues are at the heart of any explicit calculation and have to be 
dealt with. Different sets of rules can and do lead to different results, and the important 
point to emphasize is that when one 'goes' explicit, so to speak, one encounters problems. 

There is no doubt that any explicit calculation should be handled with care, and in any 
case, should be in accordance with general rules and checks j^. But in our specific situation, 
relations ()5.12j) and ()5.13|) are not of any direct help. Indeed, since we have an explicit on- 
shell expressions for Im*!!; and Im*nf to order p^, from j5J| and ()4.42j) respectively, we 
could hope to obtain a constraint from ()5.13|) that we try to check whether it is satisfied. 
Unfortunately this is not the case. Indeed, because *Ilau satisfies the second relation of 

n 

fl5.1()|l . it depends on three independent functions *!!;,* 11^ and such that j48l |: 

i 

*n,, = {s^^ - py) *Ut + py + 2 *n;) . (5.15) 

What happens is that constraint (j5.13j) . when applied to lm*Il^i,, implies only the condition: 

Im*n;(P) = 0. (5.16) 

Recall that the choice of the strict Coulomb gauge was precisely made to avoid having to 
deal with *]![ which is hence of no direct physical significance, as much is therefore 
(|5.16|) . It would still be of academic interest though to check whether indeed an explicit 
calculation, through similar explicit steps as in this work, would satisfy ()5.16|) . This would 
take us way beyond the scope of the present article. 
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C. Other explicit calculations 



Besides, it is important to emphasize that there are other exphcit estimates of the gluon 
damping rates where approximation schemes have been used, away from formal manipula- 



tions. One such estimate is that of j5j| where result (jl.3|) is derived. We have discussed in 
detail in 56[ the difference between this result and the one we carry here: ()1.3|) is obtained 
by integrating ultrasoft loop momenta < k < fi while keeping p just soft, and the infrared 
behavior is regulated via the introduction by hand of a magnetic mass directly in the static 
magnetic propagator. Result ( II. 3|) cannot be carried as it is to the limit p since it 
would give zero whereas the dampings there are finite. 



Interesting also is the result (jl.4j) from 



which is the same as (|1.3I) except that 



the constant a is determined to be equal to Nc/^ii. It is argued in [5^ that when the 
external momentum p of the quasi-particle is of order gT, the dominant term in the next-to- 
leading order self-energies comes from loop integrals involving the dressed transverse-gluon 
propagator * A^, and it is argued that such dominant terms come proportional to the function: 

Sdp) = TY. / T^*'^* W *A.(P-ir)U_^(,)^o. , (5-17) 

fco ^ ^ > 

where i indicates any of the quasi-particles. Using the structure of the spectral function 
corresponding to *At(i^), it is argued that the dominant contribution to the Matsubara 
sum in fl5.17|) is the 'static' contribution = 0. The integration over the spatial loop- 
momentum is then done for this static contribution with the magnetic mass /i taken as the 
infrared cutoff, i.e., j^"^ dk dk. The logarithmic infrared sensitivity In {nig/ n) is 

obtained as a consequence of the fact that the calculation is performed on mass-shell, and 
result ()1.4j) is said to be 'mass-shell' infrared sensitive. 

But it is stressed that the leading infrared sensitive result ()1.4j) should not be carried 
from the soft region p ~ gT where it is determined to the ultrasoft region p < ~ g'^T 
because there, the logarithmic sensitivity of Si (P) is softened and the result becomes finite. 
It is also argued that in the ultrasoft region, formally higher-order terms in HTL-summed 
perturbation may contribute to the self-energies with the same magnitude as the next-to- 
leading-order one we are concerned with in this work. 

It is clear that here too our result ()4.43|) is not in any way in contradiction with ()1.4|) 
or in fact with any other result available in the literature: we carry the calculation directly 
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in the ultrasoft region and sum soft loop momenta in the context of full one-loop-order 
HTL-dressed perturbation. 



D. Momentum expansion in scalar QED 

But we may still question ()4.43p since it depends strongly on the manner we have carried 
out the calculation. A crucial step in this latter is to perform a momentum expansion 
before the sum over Matsubara frequencies and analytic continuation to real energies are 
done. Given the fact that the final result is significantly infrared sensitive, it is legitimate to 
question the validity of this method. It may also be true that our calculations do not provide 
the equality between the damping rates of non-moving longitudinal and transverse gluons 
because intricate intermediary steps, dealing with potential infinite terms, are not handled 
with sufficient care 0]. We have put forward arguments in previous works regarding this 
issue but probably the best way to show that the method used has little to do with 

the infrared sensitivity of the damping rate is to evaluate some physical quantity with no 
momentum expansion, then perform the expansion and compare the result with that where 
the expansion is done at an early stage, before Matsubara sum and analytic continuation 
to real energies are done. It would be most suitable to make this comparison directly in 
massless QCD, but unfortunately the calculations are forbiddingly difficult. However, the 
structure of scalar quantum electrodynamics (scalar QED) is manageable enough that it is 
able to offer such an opportunity. 

The HTL-summed photon self-energies are evaluated for all u and p in hot scalar QED 
without recourse to any momentum expansion [g^. Let us focus on the longitudinal case 
and quickly review the corresponding results, with a slight editing in the notation from . 
Three regions in u and p are to be distinguished: u"^ < p^ (region I), < cj^ < Ami + 
(region II) and 4mg + p^ < u"^ (region III), where irtf, = eT/2 is the scalar thermal mass 
and e the coupling constant. The longitudinal next-to- leading order HTL-summed photon 
self-energy is found to have the following expression: 



Stt p^ 



Anie + 2ie — i — In 



uj f 2me — i{e + p) 



(5.18) 



p \2me — i{e — p) 

where the prefix 6 indicates here that the known leading hard-thermal-loop contribution has 
been subtracted for ultraviolet convergence, and e = in regions I and HI and e = 
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in region II, with 0^ = uJ^ (cu^ — — Am^) / {uj"^ — p^). The longitudinal photon dispersion 
relation is: 

n^=p^ + Ui{n,p), (5.19) 

and, up to next-to-leading order, writes: 

tf(p)-^?(p)+ , 7'i"'f' y (6.20) 
where SUi {uj,p) is the hard thermal loop given by (621 ]: 

^n,(..p)^3,„j(l-^)(l-|,„^). (5,21) 

with TJip = eT/3 the photon thermal mass, and Ui (p) the on-shell longitudinal photon energy, 
solution to ()5.19|) to lowest order, where only the hard thermal loop is kept in the self-energy. 
Using (j5.2H) . we can rewrite (|5.2(jp as: 

fi2 ( ) =2^ (uoi^p) + ^ . [ ^S*Ui {uup) . (5.22) 

3m^ + p^ — cuf 

Remember that all these results involve no expansion in p and the Matsubara sum and 
analytic continuation to real energies are done. Now we expand. The region of interest to 
us is region II, where we are allowed to perform the following expansion: 

cui {p) = l + - ^p" + O (p^) ; p = p/nip. (5.23) 



Using this, we perform the expansion of (|5.22|) for soft p. We find: 



1 - 0.368 e + ...)+( ^ - 0.0536 e + . . . ]f + 0{p'') 



(5.24) 



Here we have together the leading and next to leading orders in the coupling e. What 
remains to do is to perform the expansion in powers of before the Matsubara sum and 
analytic continuation to real energies are done, in a way similar to what we have done in 
this work for the transverse gluon damping rate. This we do and we obtain the same result 
fj5.24|) . This clearly indicates that the calculation method we use is likely to have little to 
do with the infrared sensitivity of the results. 

E. Final comments 

Trying to determine the damping rates for quasi-particles is essentially a mean of trying 
to understand the (analytic) structure of QCD in the infrared at finite temperature. Already 
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a appears at order p^. If by some mean we are able to manage an expression for the 
coefficient of p^, should we expect a This is not certain. Indeed, if we look at the 

damping rates 7±(p) for quarks in an expansion in powers of p, the zeroth-order and first- 
order coefficients are free from infrared divergences; only the second-order coefficient, that 



of p , is infrared sensitive, and only logarithmically 

7±(j9) = ^ [5.705672 T 1.056796 p+ (8.553573 - 5.968072 In /i) + O (p^)l , (5.25) 
IbiT 

where p = p/mj and jl = fi/mj with mj = a/C/ /8 gT the quark thermal mass. Recall that 
Cf = {N^ — 1) /2Nc and here, Nc = Nf = 3. The same behavior is obtained for other values 
of Nc and Nf. We are not able to explain why there is such a difference in the infrared 
sensitivity between the quarks and the gluons. 

Once the work on the dampings is finished, it would be interesting to look at uj^^\p), the 
correction to the energy at next-to-leading order in HTL-summed perturbation. One would 
have to deal with the real part of the HTL-dressed one-loop-order self-energy, somet hing 
much more difficult than dealing with the imaginary part. There is already the work (6ll |. 
where the next-to-leading order correction to rrig is estimated in the pure gluonic case using 
the longitudinal one-loop HTL-summed self-energy. The result is: 

u;l = mlil-1.8Ncg + ...), (5.26) 

and is gauge independent within the class of covariant gauges. The quantity Uq is the long- 
wavelength limit of the frequency spectrum and it is this result we have used in the text 
after (|5.7|) . As one can see, the result is finite and the HTL-summation scheme poses no 
problem. But how about trying to calculate the energy for moving ultrasoft quasi-particles? 
Is there going to be infrared sensitivity as for the damping rates? 



APPENDIX A: EXPANSION IN POWERS OF MOMENTUM 

In the main text in section four, small-fc expansions have been frequently used. They 
relied partly on the expansion of the gluon on-shell energies, the residue and cut functions 
together with their first and second derivatives. We give in this appendix the first terms of 
these expansions. We note that in actual calculations, it very often occurs that the actual 
expansion has to be pushed quite further, especially for the residue functions and their 
derivatives. We do not show these extra terms here to avoid saturation of the display. 
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Small-/c gluon energies are given by: 

)- + U ~ ^ + 6000 ^ 43120000 ~ 509600000 ~ 543312000000 ^ 
_ %k^ 88k^ _ 91617A;« 706662A;i° _ 264087844A;i2 

^ ~ ^ ^ ~ "35" ^ "375" ~ 336875 ^ 1990625 530578125 ^ ' 
Their first derivatives are given by the following expansions: 

,,,, 3k 3k^ k^ 489A;^ 79k^ 22129^;^^ 



5 70 1000 5390000 50960000 45276000000 
2677 A;^^ 

+ ...; 



381180800000 

/ (i,) _ 6A; _ 36A;=^ 176A;5 _ 732936fc^ 1413324^;^ _ 1056351376A;ii 
'^^^ ^ ~ T ~ ^5~ ^ ~l25 336875 ^ 398125 176859375 

3815200504A;i3 

^ 372246875 ^ ' 
and their second derivatives by the following expansions: 

_ 3 _ 9A;2 k^ mk^ _ 7nk^ _ 22129A;i° _ 2677A;i2 

^~5~7(r^200^ 770000 ~ 50960000 ~ 4116000000 ~ 29321600000 ^ 

6 lOSP 176A;^ 732936A;6 12719916A;8 1056351376fc^° 

uj"(k) = 1 \ 

*^ ^ 5 35 25 48125 398125 16 078125 

3815200504A;i2 

^ 28634375 ^ "" 
The residue functions are expanded as follows: 



3^ 9A;^ k^ A89k^ 711A:^° 22129A;^^ 
U ^ 280 ~ 1200 ~ 6160000 ^ 509600000 ^ 49392000000 ^ 



1_2A;2 19A;4_314A;6 736A;«_ 13717342A;i° 697037162A;i2 

' ~ 2~ ~^ ^5 375~^ ^39 5971875 ^ 176859375 

10718685944A;^^ 



+ 



1563436875 
Their first derivatives are given by: 



^ _ 9^ 1467A;^ _ 711A:^° _ 22129A;^^ 

~ 280 ^ 600 ^ 6160000 ~ 127400000 ~ 9878400000 ^ " 



_ _Ak 76k^ _ 628k^ 5888k'^ _ 27434684A;9 2788148648A;ii 

~ ~y ^ ^5 125~ ^ 539 1194375 ^ 58953125 

21437371888A;i3 

223348125 ^ "' ' 
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and their second derivatives by: 
3 9 



711k^ 



22129A;« 



8031A;i° 



J - —j^ ~ 280 ^ 200 ^ 1232000 ~ 18200000 ~ 1097600000 ~ 18659200000 

_ _4 228A;2 _ 628A;^ 5888A;^ _ 82304052A;« 2788148648^;^° 
^ ~ ~5 ^ ^5 25~ ^ ^7 398125 ^ 5359375 



+ 



21437371888A;^2 
17180625 



+ .. 



(A6) 



The longitudinal cut function expands as follows: 



1 

20' 



-k- 



1097 9tt^ 
+ 



42000 



1065078779 



8000^ ^ + 

3029563^2 



118511 68077r2 8l7r^ 

+ + „ 



12600000 5600000 3200000 

92907^^ 729a« 



k' 



+ 



+ 



363825000000 3920000000 2240000000 1280000000 
16 801915872097 8591696385777^ 594281977r^ 



k' 



+ 



19864845000000000 226380000000000 
160 1 377r6 656l7rS 



+ 



1568000000000 



+ 



128000000000 512000000000 



A(l 



^t, k) = —k 



197 



+ 



5250 
120421003 



^1 k^- 
lOOOj 

146 8 637r2 



451 



185 77r2 



28125 175000 100000 
281077r4 7297r*^ 



3789843750 30625000 17500000 10000000 

2454132521881 8570581777^ 5805097r^ 

+ 



k' 



77597050781250 110537109375 
507877r^ 656l7r« , . 



382812500 



250000000 1000000000^ 
In actual calculations, we had to expand the above cut functions till 0{k^^) 



(A7) 
The same is 
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true for the following transverse cut functions: 



(3t{l-u;uk) 



40 



176 



516128 



3157r2 36457r4 



k- 



33298514048 827024 



+ 



24721 



73811257r6 30 3 757r4 9922507r 



+ 



2148286932320768 1508396467072 3742632736 



149467781 257r8 



1025393 



12916968757r6 167442187571^ 19100812507r2 



138598879725606668288 24468480689779456 1394874080310776 
302672257031257rio 523137234375^ ^ 10172112890625^ 

38108673112 52966144687 



6446524218757r4 6257426175000007r2 



20 



97r2A; 



+ 



173056 



184 



36 457r4 1057r2 



k- 



7487094784 3381248 



255916 



73811257r6 708 7571^ 4961257r2 



+ 



323921668734976 533311520768 



6589211648 



153921104 



14946778 1257r8 



43 5 6 5 6 257re 



16744218757r4 3183468757r2 



14014147076150001664 



30 2 6 722 5 703 1 257rio 1743790781257r8 
1529299484672 11073473445148 " 



5431917214171136 1167210299260928 

+ 



10 1 721128 9 6 257r6 



+ 



2148841406257r4 195544567968757r2 



k' 



The expansion of the first derivatives of the longitudinal cut functions writes: 



dMk,uj)\ 



UJ=\—U)l 



1 

20 



2459 277r2 



42000 
3467595349 



+ 



8000^ 
107189597r2 



20027 



+ 



1837r2 



8l7r^ 



360000 1120000 640000 
1764997r^ 510371^ 



+ 



121275000000 3920000000 2240000000 1280000000 
22977822170892 7 20397448167971^ 827883277r^ 



2039744816797r2 
+ ^..^ + 



19864845000000000 75460000000000 
62 71837r6 5904971^ 



128000000000 512000000nn000 



dkl3i{k,u)\ 



ljJ = l—UJt 



1 

To 



125 9 277r2 



5250 
53836759 



k^ + 



+ 



1000^ 

18110597r2 



7937 



1568000000000 

k' + ...- 

6337r2 



8l7r^ 



39375 35000 20000 
225997r'* 51037r* 



k' 



1624218750 30625000 17500000 10000000 
4091313116071 180415270l7r2 29708197r^ 



77597050781250 36845703125 
15333371^ 590497rS 



382812500 



250000000 1000000000 



(A8) 



(A9) 



It is important to note that these expressions are not obtained by simply deriving ()A7|1 . 
One has first to calculate dk(3i{k,uj), then replace u hj 1 — ui^t, and only then perform the 
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expansion. The derivative of the transverse cut function is given by: 



-40 



2735072 



76 



319613296256 24183296 



89813 



73811257r6 



2 126 257r4 99 2 2 507r2 



+ 



5971222733143552 2602609504448 1874876804 



+ 



+ 



3814589 



29893556257r8 2583393757r6 3348843757r4 7640325007r2 
2522063079101393782784 45494694191332864 



302672257031257rio 



747338906257r8 



1031369244173616^ 



98656172204 



145623897509 



101721128906257r6 64465242 18757r4 62574261 75000007r2 



dkPt{k,uj) 



\oj=l—u>t 



-20 492544 
+ 



16 



+ 



97r2A;2 36457r4 
423589874499584 



1057r2 



35131752448 3223552 



487252 



73811257r6 



236257r4 4961257r2 



565829435392 5720261632 



45637336 



29893556257r8 86113 12577^ 33488437577^ 6366937577^ 
1 17503233176950013952 411 16968270430208 



+ 



302672257031257rio 17437907812577^ 
8361475700424704 482750388224 15508324754764 



(AlO) 



10172112890625776 2387601562577^ 1955445679687577^ 



r + .... 



The derivative over cu of the longitudinal cut function is expanded as: 



d^Pi{k,uj) 



\u!=l—UJl 



-1 

'6k 



19 977^ 



900 
245274179 



k + 



4063 



50777^ 



2777^ 



252000 56000 



+ 



^3 

64000^ 
I7OI776 

19845000000 ~ 392000000 ~ 112000000 " 128000000 
371842066997 97384005177^ 6868109177^ 



1487147772 6552977^ 



65488500000000 1078000000000 156800000000 
16629377^ 1968377^ 



6400000000 



(jJ = l—UJt 



-1 

6fc 



56 

225 



51200000000 

9772 



11507054 



k + 

200 

404403772 



2777^ 



1709 33772 ^ 
7875 ^ 3500 ~ 4000 
407777^ 1701 



310078125 6125000 875000 2000000 
26289763331 764943847772 53652377^ 



+ 



511628906250 12632812500 76562500 
3264377^ 1968377^ 



+ 



25000000 200000000 



k^ + ...; 



(All) 
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The corresponding results for the transverse cut function are: 



\u)=\—U}l 



+ 



-400 

95188 
"2976757r2_ 
4350391 



739648 



+ 



52 



1 

k 



7297r4 1897r2 

3822935800822784 



219717754112 16061392 



k + 



44286757r6 
547835807552 



1275757r4 
400625132 



1528065007r2 



17936133757r8 516678757r6 2009306257r4 
172430102335794122752 194970192021158656 



+ 



1555758096434912 



du,Pt{k,uj) 



400 



\uj=l—uit 



277r2A;3 



+ 



53248 116 



20178150468757rio 3138823406257r8 
80143211084 2836963147 1 , 

683552 ■ 



20344225781257r6 5525592187571^ 62574261 7500007r2 



7297r4 637r2 



6335234048 2425856 



44286757r6 4252571^ 29 7 6 757r2 



+ 



49834102882304 26537099264 6961211392 



14887772 



17936133757r8 172226257r<5 20093062577^ 127338757r2_ 
1078011313550000128 576925575675904 695975917125632 



20178150468757rio 149467781257r8 203442257812577^ 
1435899206656 5003479782988 " 



1289304843757r4 391089135937577^ 



r + .... 



(A12) 



We turn to the second derivatives. We have for the longitudinal cut function: 



dMk,^)\ 



1 



UJ = l—LOl 



lOA; 



+ 



73 1 277r2 



7000 



4333821991 



k + 

2000^ 
40088437r2 



144943 21577r2 



2437r^ 



1260000 280000 320000 
250297r^ 510377^ 



+ 



+ 



25987500000 980000000 1120000000 160000000 

761041 1777^ 



k^ 
k' 



8237191315679 2839378202377^ 

+ 



73573500000000 3773000000000 156800000000 



1 



5321777^ 
3200000000 
681 



+ 



590497r« 
51200000000 



277r2 



875 250 



17678 197777^ 2437r^ 

+ 



39375 8750 



5151227623 103344371^ 2315797r^ 



10000 

51037r6 



5684765625 



7656250 



^750000 1250000 



+ 



979623985177 1 7149032487r2 23548597r^ 



+ 



2586 568359375 7369140625 
51037r6 590497rS 



38281250 



+ 



6250000 100000000 



k' + .... 



(A13) 
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The corresponding results for the transverse cut function are: 



ijj=i—ijji 



4841152 



+ 



12157r4 57r2 
'9324383097777152 



1034972897792 41030672 353 

+ 



24603757r6 
8040766589504 



+ 



708757r4 
4285921012 



8757r2_ 
8089 



3321506257r8 86113 1257r6 1116281 257r4 1575007r2 

154471 11 103341084491 776 479468923614041344 



+ 



100890752343757rio 581263593757r8 
31407182749531952 399195734764 463459 



33907042968757r6 2148841406257r4 33 6 8 750 007r2 



u>=l—u>t 



763904 



+ 



12157r4 57r2 
433173355823104 



83111182336 60286976 2424 

+ 



24 6 3 757r6 
1589288370176 



+ 



708757r4 
4023629824 



8757r2 
27848 



k 



+ 



3321506257r8 2870437577^ ' 111628 12577^ ' 3937577^ 
438501832773261590528 20721255879540736 



1008907523437577^0 



830376562577^ 



+ ■ 



41103339808096256 16084864550912 3965392 



3390704296875776 21488414062577^ 6015625772 



k^+.... (A14) 



We also need the second derivatives with respect to u. We have for the longitudinal cut 
function: 



dlPi{u;,k))l 





"3 


377^1 


1 


' 51 


7I772 


977^ 




5 


40 


k ~ 


175 


5600 


1600 



k + 



395161 1715177^ 



1532777'^ 
2800000 " 
27402377^ 



I7OI776 1 ^3 
640000oJ 
656177^ 



2221367965177^ 



5250000 
413759777^ 



700000 



k' 



970200000000 1960000000 

192681536863 410415001377^ 



640000000 640000000_ 
404542729377^ 387150377^ 



36786750000000 339026688000 
11591177^ 7217177^° 



+ 



22638000000000 11200000000 5120000000 204800000000 



U)=l—U!t 

I7OI776 

+ 



377^ 



20 



k 



+ 



162 



319772 977^' 

-i 



175 700 200 



k- 



69704 255477^ 



k'+...; 
492377^ 



+ 



+ 



200000_ 

656177^ " 
5000000^ 
222639377'5 



r + 



1256517454 925954283772 

+ 



46875 9375 87500 
202860377^ 5427776 



k' 



1894921875 
696860967J 



2526562500 



15312500 5000000 



287396484375 19159765625 
218777^ 7217177^° ,7 



11304813868772 159177034377^ 
+ 



44214843750 



87500000 2000000 400000000 



(A15) 
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For the transverse cut function we have: 



748658115584 



-8000 



uji 817^2^5 



2000 27904000' 



1897r2 



+ 



1 
k 



+ 



8857357r6 
48565749544951808 

10761680257r8 
37302948014309056 



14773972 



21877r4 

5463228688 



1 



178084 188938432 



59 5 357r2 

6853501867264 



7654577^ 



229209757r2 172226257r4 310007257r6 
278379176783 3300542659852 

k 



12206535468757r6 



46930696312507r2 7735829062571^ 

309677336393249792 26055 1 1 136700461694976 

+ ■ 



dl(5t{Kuj)\ 



■1000 



209254893757r8 



1000 1088000 



+ 



+ 



U)=\—U!t 



13866631168 



8l7r2A;5 



637r2 21877r4 



1 

fc3 



12106890281257rio 
233368 29452288 



1 



8857357r6 
476299021778944 



41076832 86458213888 
+ 



595357r2 255 1 57r4 



456797978624 



10761680257r8 
43774954685857792 
12206535468757r6 



76403257r2 

k + 



12055837571^ 103335757r6 



7827935167072 4077064460288 

+ 



+ 



23465348156257r2 
3273422343766016 
23250543757r8 



257860968757r4 
14262918917738463232 
12106890281257rio 



(A16) 



APPENDIX B: THE OTHER TERMS 



In the main text, we have shown how to determine the Zgll contribution to the coefficient 
of order in 7^ (p). Here we give the expressions of the other contributions. Their derivation 
follows the same steps as in the main text. Whenever there is an infrared divergence, it 
is extracted as explained there. All hght-cone divergences are brought under control as 
explained too. Recall that all is in unit of irig and we use the same notation as in the main 
text, except that fi — /i/mg is directly used instead of /i. 

We have: 

/ V— [4 + 93A;^ + ISk'^ + (8 + 78k^)cu - (140 + 78k'^ + 9Qk^)u;^ + (48 - 180k'^)u;^ 

/n. .^.,2N 4 ,^^5 61 / 0.7505272862 6.404458463 
- (31 - 126A;^)a;^ + 102a;^ - 540;"^] ptp'i = h 



+ 97.06037391n/i + 73.01499652. (Bl) 

The infrared behavior here is,l/jx^ together with and In Jx. Because of this, we have good 
convergence starting from I — 0.2 but the convergence is lost from I — 0.025 downwards. 
This is not a problem because I does not have to be so small. The corresponding double 
integral is finite. We make the change of variable k — 1/k. To get the number, we have 
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first to push the boundary lo — 1/k to cu — 1/k — e and the hmit k — 2ioK — 2 — e'with 
the condition e' > 2e (in order for the arguments of the logarithms in the cut functions 
to remain positive). With this, there is good convergence to the number —4.96229365... 

starting already at e = 10~^. 
The next term is: 
1 



I V— -18A;2 + 122A;^ - (144 - nAk'^)ujuo' + ( - 138 
J ujuj' [ 



2 12 

UJ UJ 



174 



UJ^UJ'^ 



^0 



PtPt = 



37r2 



In /i - 9.45862819 . (B2) 



Here I — 0.2,0.1 is sufficient. There is in fact good convergence: we have a good num- 
ber already for £ — 0.3. We note here that we have to replace in dxk'^^tPtf the upper 
boundary 1 by 0.9999... because the integral does not converge numerically in the ultraviolet 
(analytically it is convergent). However, for fixed £, say i = 0.1, there is stability of the 
integral when adding nines. The corresponding double integral does not pose a problem. 
Next we have: 



V- 



9 3 \ 9 

2^+4* +r+ 



6 



f 57 



57 



_3_ 33\ 2 _ o_ , 



9 u 
+ ttC^ + 



27 



tUJ 



Gp; = ^In/i- 4.904610308. 



(B3) 



k^ Ak^ 

For the double integral, one has to regularize at a; = 1 — k: one shifts the boundary to 
UJ — 1 — k — e. We have good convergence starting at e = 0.0001. 
The next contribution is: 



pi 



9 o 15 2 6 3 39 4 9 5 
h 3cc; - —UJ^ + —UJ^ + —rUJ^ - —^uj^ 

2 A;2 fc2 2k^ 



0(9fcp; = 4.464952656. 



(B4) 



This contribution is finite. The corresponding double integral has initially an apparent 
light-cone divergence which is handled as explained in the text. 
The next term is: 



/ 



V— [mk + 14A;^ - ?>k^ - {hAk + l&k^ - 6k^) u + {2k + lAk^) + {Sk - 12k^) 

,,,4 ,,51 ^ , 1000 /34120 54701440\ 1 
-life' + 6k.=] p,a,A = + (y89^ + ^ 

/4277584 1938884128 1278453185024\ 



UJ 



+ 



+ 



+ 



V595357r2 7654577^ 88573577^ J 

Here too there is a l//i^. But there is no fight-cone divergence. 



In /x+ 1477.282103. (B5) 
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The next term is: 

/ '^uu? ^""^^^ ~ ''^^^ ~ ^^^^ ^ "^^^^ ~ ^^^^^ ^ ~ ^^^^ ~ "^^^^^^ ^ ^^^^ ~ ^^^^-^ 



400 



/ 12136 8r)6G912\ 1 
+ 



317576216 2290058432 1411647241216\ 
5+ o + .SI J^ . . + 



37r2/i4 \^ 637r2 7297r4 J 

ln/x + 1109.862748. 



(B6) 



99 2 257r2 ' 12757571^ ' 14762257r6 J 

The stabihty band here is a httle narrow: I — 0.1,0.045. Also, the double integral has an 
apparent hght-cone divergence and is treated as usual. We note that in order to have a 

convergence in the double integral, fcmax = 400. 
The next term to discuss is: 



V- 



UJUJ' 



32k^ + 20A;^ + (112A; + 124A;^) cu + - 36k - 116/c^^ cu^ - 

PtdkP't ■- 



108 \ o 

— + 232A;J uu^ 



12 



108 



+ ( y + 156A; ^ 



There are no hght-cone divergences here. 
The next integral is: 



60 



k 



928 



277r2 



In/x- 1.97446447. 



(B7) 



y"p^ (^-3k+^u^^ Qdlp't = -3.245787103. 



(B8) 



No infrared and no light-cone divergences. 
The next contribution is: 



/ V—^ [30k^ + 4k^ - 2A;^ - (24^;" + 8A;^) uj - {Ak^ - 4A;^) cu^ + 8^00^ - 2kW] pidlp'^ 

ln/i+6810.290071. (B9) 



19364608 16 



7297r4 37r2 



176 6906341632 8279783182336 
4 



257r2 1275757r4 



14762257r6 



The double integral has a light-cone divergence that is treated as usual. 
The next contribution is: 



/ 



V—^ W + 16A;6 + \6k^uj + (8A;2 - 32A;^) + 16A;2^^1 p^dlfJ^ = 
ujuj' '- \ / J 97r^//* 



123376 170390272 

+ 



1897r2 



21877r4 



p^ 



515581456 11866583750656 

+ o + 



992257r2 



14762257r6 



+ 



4181643904' 



425257r4 



In /i - 8137.840681 



(BIO) 
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Because of the presence of l//i^ and because we are requiring more precision, the stabihty 
band for the single integral is narrow: £ = 0.1,0.05. The double integral is light-cone 
divergent and handled as usual. 
The next contribution is: 

/ V—^ [8k'^ + (20fc2 + 20k^) u + {l2 + 12k^) cu^ - (l2 + 32P) - 12uj^ + 12tu^] ptdlp[ 
J '^'^ 

4400 , 

In/i- 10.70561628. (Bll) 



277r2 

Generally, the transverse terms are 'smoother'. The double integral has no light-cone diver- 
gence. 

The next contribution is: 

/•+00 /"+00 /•+00 _Q /I 2 _ ^2\2 

/ dk \ diu duj' ^— ^ — ^Gp; 9^5(1 -a; -a;') = 4.65989517. (B12) 

Jp, J-oo J ~oo k UJ 

No infrared or light-cone divergences, though this term is highly 'uncommon'. 
The next contribution is: 

p+oo p+oo p+oo KAh'i „ 8000 

/ dk duj duj'- -cj2)'p,p;a^5(l -cj-cu') = In /i- 268.9616004. 

Jj, J-oo J-oc Stt^ 

(B13) 

The double integral has a light-cone divergence which is treated as usual. 
The next contribution is: 

p + OO p + OO p + OO 1 QO 1 fiDO 

/ dk duj duo' [k^ + uouo') ptp[dj il - uo - uo') = ^In/i- 34.46857932. 

J-oo J-oo 37r2 

(B14) 

The double integral has no light-con divergence and converges fairly well. 
The next contributions are: 

r (^2 _ J^2^ ^ ^_g. f iuj^ - k^) pt = 3.0. (B15) 

J OJUj' ^ ' J Uj' ^ ' 

[ V^e{uj)pid^ [5{k~u:)~5{k + uj)] = -3.0. (B16) 
J ^' 

There are no light-cone problems with these integrals. Note that the values here are 'exact' 
though the integration is done numerically. We have no analytical check of this. 

These are all the terms that intervene in the determination of the one-loop HTL-dressed 
transverse-gluon self-energy. When added to the one discussed with more detail in the main 
text, we obtain the result displayed in ()4.42p . 
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